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ABSTRACT 


v^n 3 - 3 - 1, 

A body in the upper atmosphere or in space will acquire an electric 
charge, or potential, which must be known to determine the motion of 
micrometeoiites, the dragon earth satellites, and to assess the behavior 
of certain experiments on satellites. 

The equation for ion and electron currents to a sphere are availa- 
ble in the literature for small bodies. For large bodies, an estimate of 
the influence of the plasma sheath is required in an attractive field. 
Poisson's equation has been solved numerically for high-velocity 
spheres, and the ion current obtained by an analysis of the ion's dis- 
tance of closest approach. Photoemission is an important charging 
mechanism for bodies in sunlight. Measurements of photoelectric yields 
'are reviewed and compared with photocurrents measured above the at- 
mosphere. Secondary electron emission upon energetic particle impact 
may also be an important mechanism, especial in the earth's radiation 
belts. The effects of cosmic rays, radioactivity, thermionic and field 
emission, collisions with dust grains, and the influence of radio- 
frequency electric fields are generally negligible. A magnetic field 
induces a potential gradient in a moving body; in addition, the restric- 
tion of electrons to a spiralling motion along the field line decreases 
the body's effective collection area. 

These charging mechanisms are evaluated for conditions in the 
upper atmosphere and in interplanetary space, and are combined into 
expressions from which the equilibrium potential may be determined. 
In the ionosphere the equilibrium potential is typically a few tenths of 
a volt negative. At higher altitudes the potential may become positive 
in the sunlight as photoemission predominates over positive ion col- 
lection. In the earth's magnetosphere the potential is sensitive to the 
ratio of electron flux to photoemission, and may vary widely. Positive 
values are limited to a few volts, but large negative values are possible. 
In interplanetary space positive potentials due to both photoemission 
and the solar wind protons are expected. 

The equilibrium potential of the satellite Explorer VIII has been 
measured in both darkness and sunlight. There is general agreement 
with theoretical values at higher altitudes. At low altitudes the meas- 
ured potentials are more negative than anticipated. This is shown to be 
due to a radio-frequency plasma impedance experiment carried on the 
satellite. 

Oju^ 0 ' 


iii 



LANK PAGE 



ACKNOWLEDGMENTS 


I would like to thank the National Aeronautics and Space Administra- 
tion for their support of my graduate work, and for their permission to 
use data obtained while at the Goddard Space Flight Center for the thesis. 
I would especially like to thank Dr. John F. Clark for his advice and 
guidance during the preparation of the thesis, and also my colleague Dr. 
Eugene J. Maier for his careful reading of the manuscript and his help- 
ful comments. I also would like to acknowledge the helpful discussions 
with my other colleagues: Mr. Robert E. Bourdeau, Mr. John L. Donley, 
Mr. Gideon P. Serbu and Mr. Ballard E. Troy. 

I am indebted to Mr. David Quinn for his many hours of data plotting, 
calculating, and computer programming. Mr. Daniel Elliot also assisted 
in some computer programming and Mr. William Zadell in some of the 
plotting and calculations. The staff of the Physical Science Laboratory 
of the New Mexico State University assisted with the Explorer VIII data 
reduction. 

Finally, I would like to thank my wife, June, for her encouragement 
and for her help in proof-reading the final text. 


v 



LANK PAGE 



TABLE OF CONTENTS 


Page 

ABSTRACT iii 

ACKNOWLEDGMENTS v 

LIST OF TABLES ix 

LIST OF FIGURES xi 

Chapter 

I. INTRODUCTION 1 

1. The Concept of Equilibrium Char t i 1 

2. Applications of Equilibrium Charge 

Determination 6 

II. HISTORICAL SURVEY AND PRESENT STATUS OF 

THE PROBLEM 13 

1 . Historical Survey 13 

2. Present Status of the Problem 21 

III. COLLECTION OF ELECTRONS AND IONS 2 3 

1. General Considerations 2 3 

2. Electron Collection 24 

3. Positive Ion Collection 27 

4. Effects of the Satellite Wake 33 

IV. PHOTOEMISSION AND SECONDARY EMISSION 37 

1. Photoemission 37 

2. Secondary Emission of Electrons Upon Electron 

Impact 41 

3. Secondary Emission of Electrons Upon Ion 

Impact 48 

V. OTHER CHARGING MECHANISMS 55 

1. Discharge Time for a Body in a Plasma 55 

2. Cosmic Rays 56 

3. Radioactivity. . . 57 

4. Thermionic Emission 59 

5. Field Emission 59 

6. Collisions with Dust Grains 60 

7. The Effect of Radio Frequency Fields 60 

vii 



TABLE OF CONTENTS (Continued) 


Chapter 


VI. MAGNETIC FIELD EFFECTS 

1. The Induced V X B Potential Gradient . . 

2. Effect of a Magnetic Field on the Direct Collection 

of Particles 

VII. CALCULATION OF EQUILIBRIUM POTENTIALS 

FOR VARIOUS ATMOSPHERIC ENVIRONMENTS . . . 

1. General Equations and Computational Procedure. . 

2. Expected Equilibrium Potentials in the Ionosphere 

for Large and Small Bodies 

3. Expected Equilibrium Potentials in the Earth's 

Magnetosphere and in Interplanetary Space 

VIII. EQUILIBRIUM POTENTIALS MEASURED ON THE 

EXPLORER VIII SATELLITE 

1. Description of Experiment 

2. Experimental Results 

3. Conclusions 

APPENDIX A - ION CURRENT TO A NEGATIVELY - 

CHARGED MOVING SPHERE 

APPENDIX B - SUMMARY 

REFERENCES 


Page 

63 

63 

64 

77 

77 

86 

97 

105 

105 

110 

123 

125 

135 

139 


viii< 



LIST OF TABLES 


Table Page 

I. Summary of appropriate equations 33 

II. Secondary electron yields for low-energy incident ions ... 50 



LANK PAGE 



LIST OF FIGURES 

Figure Page 

1. Phase Diagram Illustrating Typical Dependence of 

Current on Charge . 3 

2. Phase Diagram Illustrate. g the Effect of a 

Negative Resistance Mechanism 5 

3. Phase Diagram Illustrating Oscillatory Behavior 

of Satellite Charge . 7 

4. Ion Current to a Large, Negative, Moving Body 32 

5. The Wake Behind a Rapidly Moving Body in 

the Upper Atmosphere 34 

6. The Solar Spectrum and- Photoemission Yields 

from 12 2 - 10 4 A° 39 

7. Experimental Current-Voltage Curve for Photoemission 

from Tungsten by Hinteregger, et al. 42 

8. Second - , ry Electron Yieiu for Electron Impact on 

Tungsten and Aluminum 44 

9* Back-Diffusior Coefficient for Tungsten 

and Aluminum 46 

10. Total Secondary Electron Yield for Tungsten 

and Aluminum 47 

11. Energy Spectrum for Secondary Electrons ......... 49 

12. Secondary Electron Emission from Ion Impact at 

Low Kinetic Energies 51 

13. Secondary Yield of Electrons for Proton Impact at 

Energies Abov 3 1 Kev 33 

14. The Effect of a Magnetic Field on Electron Collection 
by a Cylinder; a) The Geometry of Collection, 

b) Integration Area in the r -p Plane ^7 

15. The Electron Current to a Cylinder in a 

Magnetic Field. ' u 

16. The Geometry of Electron Collection by a 

Sphere in a Magnetic Field 72 _ 



LIST OF FIGURES (Continued) 

Figure Page 

17. The Electron Current to a Sphere in a 

Magnetic Field 75 

18. Equilibrium Potentials for a Small Body at Rest 

Without Photo- or Secondary Emission 78 

19* Equilibrium Potentials for a Small Body at Rest 
as a Function of Photoemission and/or Secondary 
Emission 80 

20. Equilibrium Potentials for a Body at Rest and Without 

Photoemission as a Function of its Size 82 

21. Equilibrium Potentials for a Small Moving Body 

Without Photoemission 84 

22. Equilibrium Potentials for a Small Moving Body as a 

Function of Photoemission 85 

23. The Effect of Size on the Equilibrium Potential of a 

Moving Body 87 

‘24. Equilibrium Potentials for a Large Moving Body as a 

Function of Photoemission 88 

25. Cold Model Ionosphere 90 

26. Warm Model Ionosphere 91 

27. Equilibrium Potentials for a Small Body in Cold and 

Warm Model Ionosphere 93 

28. Equilibrium Potentials for a Large Body in Cold and 

Warm Model Ionosphere 94 

29* Typical Energetic Particle Fluxes in the Earth's 

Magnetic Equatorial Plane 96 

30. Positive Equilibrium Potentials for a Body in the 

Magnetosphere 100 

31. Positive Equilibrium Potentials for a Body in 

Interplanetary Space 102 

xii 



LIST OF FIGURES (Continued) 

Figure Page 

32. Explorer VIII Satellite, Showing Component 

Locations 106 

33. Block Diagram of the Electron Trap on Explorer VIII. 108 

34. Typical Current- Voltage Curve Obtained with the 

Electron Trap m 

35. Measured Electron Temperatures for Magnetically 
Quiet Days; Lines Indicate Mean Values Used in 

7 1 I O 

Calculations 

36. Measured Electron Densities for Magnetically Quiet 

Days; Line Indicates Values Used in Calculations ... H4 

37. Measured and Calculated Equilibrium Potentials for 

Explorer VIII on Magnetically Quiet Days 

38. Measured and Calculated Equilibrium Potentials for 
Explorer VIII on Moderately Magnetically Active 

Days 1 1 6 

39. Measured and Calculated Equilibrium Potentials for 

Explorer VIII on Magnetically Disturbed Days 117 

40. The Response of a Radio-Fi ..c. tency Probe in a 

Plasma 122 

Al. The Potential in Front of a Negatively-Charged 

Moving Sphere; Contours of the Minimum Distance of 

Closest Approach are also Shown for Various Ion 

Energies 126 

A2. Equivalent Potentials for a Given Potential 

Distribution and Energy u, but for Various Impact 
Parameters . . 129 

A3. Equivalent Potentials for Various Potential 

Distributions Denoted by (x Q » y 0 )» but for a Given 

Energy and Impact Parameter * 130 


xiii 



CHAPTER I 


INTRODUCTION 

1. The Concept of Equilibrium Charge. A body in the upper atmosphere 
or in space, such as a satellite, a meteor, or a dust grain, is continually 
bombarded by environmental particles of which a proportion will be 
electrically charged, such as electrons, ions or cosmic rays. When 
such an encounter occurs there is, in general, a transfer of charge either 
to or from the body. Ott jt processes may also occur that can effect a 
charge transfer. Incident photons of sufficient energy will induce the 
emission of photoelectrons. Under the proper conditions thermal emis- 
sion or field emission of electrons could occur. Other mechanisms of 
charge transfer include radioactivity of the material composing the body, 
secondary emission of electrons and collisions with other bodies such 
as micrometeors (dust grains). 

The rate at which charge transfer proceeds for a given mechanism 
depends both upon characteristics of the body such as its surface area 
and material and on the environmental conditions such as the number 
density of charged particles. In particular, this rate of change of charge 
on the body depends on the net charge already residing on the body. This 
is merely saying that the motion of charged particles in the vicinity of 
the body is influenced by the elecvric field arising from the charge dis- 
tributed on the body. A positive charge will attract electrons and repel 
ions. Secondary or photoelectrons may not escape but may return to the 
body. 
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These statements may be summarized in a differential equation 
for the charge on such a body. If we let I denote the total current to 
the body, we have 


I = ^ = f (0. body characteristics, environmental conditions) (1.1) 
d t 


where Q is the total charge on the body. A complete understanding of 
the charge as a fu* ction of time requires a solution of (1.1) where the 
right hand side expresses the sum of the various currents to the body 
as a function of Q and t. However, it turns out that typical charging 
or discharge times are small compared with the time in which signifi- 
cant changes in the environmental conditions usually occur. For ex- 
ample, the period of revolution or tumbling of a satellite will be on the 
order of a second or greater, whereas discharge times will be on the 
order of milliseconds or less, as will be shown later when the various 
mechanisms are discussed in detail. Consequently, it is a good approx- 
imation to assume that the environmental conditions remain constant 
during the time that it takes for a body to acquire a charge; hence, the 
right hand side of equation (1.1) will not be explicitly dependent upon 
time. One exception to this assumption occurs when the effect of radio 
frequency voltages on current collection is considered, as discussed in 
Chapter V. 

W*ith this approximation it is convenient to discuss equation (1.1) 
by referring to a "phase diagram" similar to the kind employed in de- 
scribing the behavior of oscillatory systems. Figure 1 is such a plot 
of the current I versus charge Q for a hypothetical but typical body in 
space. Such a curve is typical, for instance, of the case in which the 
positive current is due to positive ion collection and the negative cur- 
rent is due to electron collection from the environmental plasma. 




Figure 1. Phase Diagram Illustrating Typical Dependence 

of Current on Charge 


3 





The behavior of the system can be described by the motion of the 
representative point P. For positive I , P will move to the right, as 
shown by the arrow indicating an increase in Q. For a negative current 
P will move to the left. When P arrives at the point on the abscissa 
marked Q E its motion will cease; the current is zero and hence Q re- 
mains constant: Q E is the so-called equilibrium charge. It is important 
to note that this is a point of stable equilibrium, since if P is displaced 
from Q e the direction of the arrows is such as to restore the system to 
equilibrium. 

The curve in Figure 1 is also typical of most natural current mech- 
anisms in that it is monotonic - the slope of the curve is everywhere 
negative. A negative slope indicates a positive (but not necessarily 
linear) resistance between the medium and the body, with a steeper slope 
indicating a lower resistance. 

The existence of a stable equilibrium point depends upon this re- 
sistive character of the charging mechanisms involved. To see this, 
consider a hypothetical mechanism exhibiting 'negative ' resistance. 

Such a mechanism is physically feasible although unlikely. For example, 
a material with secondary emission characteristics can be imagined 
with a secondary emission yield which increases at some threshold 
energy of the incident primary electrons to a value greater than unity 
with increasing energy. If, in addition, the secondaries were all emitted 
with an energy distribution centered at an energy sufficient for escape, 
one would have a phase diagram similar to that shown in Figure 2. 

This system has three possible equilibrium points as indicated by - 
Qi , Q 2 , Q 3 . Only two of these are stable, however. The point Q 2 on the 
negative resistance portion of the curve is unstable. Even though dQ/dt 
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= dQ/ dt 




is zero at Q 2 , a slight fluctuation of the charge will cause the system 
to move to one of the two stable equilibrium points as indicated by the 
direction of the arrows. 

The question arises: might not a current mechanism be possible 
such that an oscillatory solution could be obtained? It is difficult to 
think of any natural mechanism of this type; but it would be simple to 
construct such a device which could be placed in a satellite. For ex- 
ample, an electric field meter to sense the polarity and magnitude of 
the satellite charge could be combined with an ion gun in such a way 
that the curve of current versus charge of Figure 3 would be obtained. 
The ion gun would be turned on when Q = Q 2 and turned off when Q = 

Qj. If the initial (non-equilibrium) charge were sufficiently negative 
the equilibrium point Q E would never be reached; the ion gun would 
cycle on and off as the representative point P travels around the closed 
loop ABCD. This demonstrates the fact that for an oscillatory solution 
of (1.1), when the righthand-side does not explicitly depend upon time, 
the current must be at least double-valued as a function of Q . In other 
words, the current mechanism must involve at least two states in such 
a way that the system can alternate between the two or more states. 

2. Applications of Equilibrium Char k lations. Knowledge of 

the equilibrium charge on a body is irr. o several areas of in- 

vestigation. In the ionosphere and also in interplanetary and interstellar 
space the flux of ions and electrons to a body such as a dust grain con- 
stitutes a loss mechanism for the charged particles in the medium.^ 
Electrons striking the body are usually captured, while ions in striking 
the body ordinarily pick up an electron and rebound as a neutral atom 
or molecule. Since the fluxes of ions and electrons are influenced by 
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Figure 



ION GUN OFF 
ION GUN ON 


3. Phase Diagram Illustrating Oscillat 
of Satellite Charge 
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the grain's charge, this information is necessary to determine the 
ionization balance of the medium. 

Spitzer has shown that photoemission from interstellar grains is a 
source of kinetic energy for the interstellar medium. 2 Since the amount 
of energy carried avyay from the grain depends upon the grain's charge, 
this quantity needs to be estimated. 

In interplanetary space the motion of micrometeorites, if they are 

charged, will be influenced by the interplanetary magnetic field originat- 

3 4 5 

ing from the sun,' * and also by the earth's magnetic field. 

The electric charge on a satellite increases the atmospheric drag 

by attracting ions of opposite polarity to the satellite which otherwise 

would not have impacted. Momentum is also transferred to ions which 

do not impact directly but ere deflected by the satellite's electric 

L 7 Q 

field. 0 ’*’ 0 It is also possible that the motion of a charged satellite 

through the ionosphere will excite plasma waves which can carry away 

9 10 

energy from the satellite. ’ 

The charge on a rocket or satellite may have an important influence 
on the behavior of experiments designed to measure the properties of 
charged particles in the atmosphere. Johnson and Meadows have dis- 
cussed the effect of a vehicle charge on the performance of a rocket- 
borne ion mass spectrometer.il The interpretation of data from ion 
and electron traps and from Langmuir probes must take into account 
the effects of the potential between the vehicle and the medium. 12,13 
Conversely, the vehicle potential may be determined from the charac- 
teristics of the data from such experiments, as will be illustrated 
later. 

It is frequently advantageous to discuss the charge on a body in 
terms of the corresponding potential with respect to the surrounding 
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medium. The charging currents to a body depend upon the body's poten- 
tial rather than charge; also, the potential is more easily measured on 
a satellite. The relation between the equilibrium charge and potential - 
i.e. the body's capacitance - is, in general, a strong function of the en- 
ironment, which is another way of saying that a body in a plasma behave, 
quite differently electrically from a body in a vacuum. For example, 
consider Poisson's equation in spherical coordinates for the potential 
p, where the space charge p , is given by (See also Appendix A.) 


p = N 0 e(l - e* e/kr ) 


(l.ii) 


Here N 0 is the ion or electron density at a great distance, e is the ele- 
mentary charge, k is Boltzmann's constant and T the plasma temper- 
ature. The first term gives the ion density and the second the electron 
density in front of a fast moving spherical satellite with a negative 


potential, d> a . If — 
R kT 


< < 1, we may write Poisson's equation as 


6 2 <p 2 dcp _ N o e 4>e 


dr a 


r dr 


«o kT 


(1.3) 


where e 0 is the permittivity. 

The solution of (1.3) is the well-known Debye potential 


e -(r-R)/L 

r 


(1.4) 


where R is the satellite radius and L ia the Debye length, 


L = 


'e A kT 


e 2 


. 1/2 


(1.5) 


9 



The "capacitance" may be determined from the field at the surface 


of the satellite and the potential <^ . We obtain 



c. ('.« 

This is simply the capacitance for two concentric spheres with a sep- 
aration distance L . Thus the Debye length, L, gives the screening dis- 
tance or sheath thickness about a charged body in a plasma. 

The subject of this investigation is the value of the equilibrium 
potential of a body in the uoper atmosphere and in interplanetary space 
as a function of environmental conditions and various body characteris- 
tics. Only naturally occurring mechanisms will be considered, and it is 
anticipated that for a body in a given environment there will exist 
only one stable equilibrium point. This investigation will be concerned 
mainly with simple bodies in the sense that the body surface is consid- 
ered to be a good conductor. Clearly one can imagine a satellite con- 
sisting of two metallic portions connected by an insulator such that each 
portion would reach its own independent equilibrium potential. Indeed, 
such "double probes" have been treated in the literature and flown in 
the upper atmosphere. 14,15 Recently some satellites have been coated 
with thermal blankets which are also good electric insulators, so that 
each point on the outer surface would reach its own equilibrium poten-" 
tial determined by the local conditions, orientation of the surface at 
that point, etc. 
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Chapter II surveys what has been done on this problem and indicates 
its present status. Chapters III, IV and V contain discussions of various 
mechanisms of charge acquisition, and Chapter VI examines the effects 
of a magnetic field. Chapter VII combines the results of Chapters III 
through VI to calculate expected potentials for a body in various environ- 
ments. Finally, Chapter VIII presents the results of measurements of 
satellite potential made on Explorer VIII and compares them with the 
predicted calculations. 
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CHAPTER II 


HISTORICAL SURVEY AND PRESENT 
STATUS OF THE PROBLEM 

1. Historical Survey . Apparently the first paper to discuss the problem 
of the equilibrium electric charge to be expected on a body in space was 
one by Jung in 1937, entitled "The Origin of Solid Particles in Inter- 
stellar Space. Jung obtained equations for fluxes of positive ions or 
electrons to an interstellar grain as a function of the potential of the 
grain. He concluded that in interstellar space the effective processes 
are photoemission and electron accretion. An efficiency (yield) of 100% 
was assumed for photons with sufficient energy to remove an electron 
completely from the solid particle. He arrived at values for the poten- 
tial of about one to a few volts positive. 

Spitzer in 1941 took up the problem in a paper, "The Dynamics of 

17 

the Interstellar Medium." He showed that the effects of cosmic rays 
on charging an interstellar dust grain are negligible. The effects of 
grain collisions with protons and electrons along with photoemission 
were considered. Using an efficiency of 10“ , he concluded that photo- 
emission does not predominate in its effect on the grain charge in inter- 
stellar space. In H II regions (where hydrogen is ionized and n e equals 
about 1/cm 3 ) he obtained a grain potential of about -2 volts independent 
of its composition or radius. In H I regions where n e is smaller (10“ 3 / 
cm 3 and hydrogen is not ionized) the effect of photoemission is notice- 
able and reduces the potential by about 50% but does not make it positive. 
The charge on a grain in the vicinity of a star was not considered. The 
effects of dust in "de -ionizing," (providing a recombination surface for 
ions) was mentioned but not discussed in detail. 
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Cernuschi criticized Jung for using too high a value for the photo- 
electric yield, and both Spitzer and Jung for assuming metallic grains 

1 8 

only. Cernuschi claimed that not every electron incident on the par- 
ticle surface is captured, especially for dielectric substances, and 
assumed a sticking probability of 10% for small negative potentials. In 
H II regions the grains are then slightly positive near stars and slightly 
negative in regions far from stars (-0.7 volts). In H I regions the poten- 
tial is positive — about 0.8 volts. For dielectric materials the potential 
is very nearly zero or perhaps slightly positive. 

Spitzer treated the subject very thoroughly in "The Temperature of 

2 

Interstellar Matter," considering the following factors in detail : 

(1) The sticking probability for electrons, which he takes to be between 
0.1 and 0.5, depending on the nature of the substance; (2) The neutral- 
ization probability for ions, which because of lack of information could 
be anywhere from 1 to 10“ 4 ;(3) The fraction of available photon energy 
converted to photoelectron kinetic energy. This involved an average 
over the photoelectron kinetic energy spectrum. He arrived at a figure 
of 0.55 for this fraction representing a value midway between the em- 
pirical value of 0.45 and the theoretical value of 0.67; (4) The threshold 
frequency which corresponds to the work function and is dependent on 
the nature of the substance; (5) The photoelectric efficiency or yield 
which is a function of both the incident frequency and the substance; 

(6) And finally, the relative absorption cross section for photons which 
differs from the geometrical cross section for particles with diameters 
comparable to the photon length. He obtained equations which can be 
used to determine the equilibrium electric charge involving all these 
factors as well as the density and frequency distribution of interstellar 
photons. 
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In a companion paper Spitzer and Savedoff in 1950 concluded that 

if the sticking probability for electrons is equal to the neutralization 

probability for ions, then in H II regions the potential V = -2.2 ( T/ 10,000°) 

19 

volts for non-metallic substances. If the neutralization probability 
is less than 10" 4 , then the potential may reach -3 volts where field 
emission of electrons is likely to occur. For highly photosensitive 
materials (which are unlikely) positive charges are expected, especially 
close to hot stars. Metallic surfaces have moderate potentials as com- 
puted by Cernuschi. In H I regions low temperatures for the gas yield 
low potentials if photoemission is weak - only a few electrons per grain. 
Especially sensitive or metallic grains may have positive potentials up 
to 10 volts if photoemission is important. 

Johnson and Meadows postulated a negative rocket potential of 20 
volts above 120 km to explain the results of an ion mass spectrometer 
experiment flown in 1954.^ They suggested that the potential could be 
due to energetic electrons or to absorption of X-rays and subsequent 
ionization of the gas evolving from the rocket. 

The first calculation of electric charge on a macroscopic body was 
apparently made by Lehnert in 1956.^ Anticipating the orbiting of earth 
satellites, Lehnert took into account the increased positive ion current 
to the satellite on its forward surface due to the high satellite to ion 
velocity ratio. The resulting potential was about -0.7 to -1.0 volts. 
Photoemission may change this value depending on the type of surface. 
The electric field caused by polarization of the satellite in the earth's 
magnetic field was said to be small compared to the field from the net 
charge on the satellite. The satellite was shielded electrically at about 
bne Debye length. 
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In 1956 Singer discussed the charge on micrometeorites near the 

5 ... 

earth. He assumed a power law for the photon energy distribution and 

balanced photoemissiori against electron accretion. He computed a 
potential of -8 volts in the dark (neglecting the ram effect on the positive 
ions). Assuming a mean energy of 1.5 volts for electrons and a photo- 
electric yield of unity, he obtained a potential of about 100 volts positive 
for particles in the sunlight. 

In 1957 Opik calculated the expected charge on interplanetary dust 
20 

particles. He assumed that the radiation effective for photoemission 
comes largely from the solar corona rather than the photosphere. 

Rather than estimating photoemission currents from assumed yields 
and work functions, he used Saha's equation for an ionization equilibrium 
applied to a diluted solar corona and the solid grains. He arrived at 
values of grain potential ranging from 50 to 220 volts positive for inter- 
planetary electron densities from one to 600 per cm 3 . 

In 1957 Jastrow and Pearse estimated the charge on a satellite in 

7 

order to find the additional drag. They neglected photoernission but 
took into account the ram effect of the satellite's velocity on the positive 
ion current plus the attraction of a negative satellite for positive ions. 
The equilibrium potential was computed to be from -10 volts on the 
night side of the earth to -60 volts on the day side because of the high 
energy assumed for the electrons (1.5 volts with a tail at higher ener- 
gies). They showed that the satellite is effectively screened at a dis- 
tance of a few Debye lengths. 

21 

In anticipation of the first Soviet Sputnik, Gringauz and Zelikman 
in 1957 discussed the distribution of charged particles around a satellite 
and derived an equation for the equilibrium satellite potential taking 
into account the satellite's velocity and photoemission: 
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0 = - 


kT 

e 


e 



( 2 . 1 ) 


Herek is Boltzmann's constant, T e the electron temperature, e the 
value of the electron charge, and I e , I + and 1^ the electron current, 
ion current and photoemission current to the satellite. The induced 
potential gradient caused by the satellite's motion in the earth's mag- 
netic field was estimated. 

In a companion paper Imyanitov discussed the problems of meas- 

22 

uring an electric field in the ionosphere from a satellite. The field 
due to the charge on the satellite must be eliminated; its magnitude 
was estimated at several volts/cm by computing the satellite potential 
in a manner similar to that of Gringauz and Zelikman and estimating 
the sheath thickness from plasma probe theory. 

Fred Whipple in 1958 used Spitzer's method of computation and 

23 

Hinteregger 'c data on sc a* ultraviolet flux and yield. With a mean 
wavelength of 1000 A and an efficiency of 0.2, he obtained a photoemis- 
sion rate of 5.7 X 10 10 electrons/cm 2 sec (9.1 X 10”® amp). If the 
electron temperature is 500,000°, then the potential is zero for meteoric 
dust if the electron density is 130/cm 3 . 

Chang and Smith in 1959 derived an expression for satellite poten- 
tial by balancing the positive ion current (simple ram expression plus 

24 . . 

a first order correction) against the electron current. Photoemission 
may also be included, but the authors concluded that its effect was neg- 
ligible. However, their expression for photoemission is incorrect in 
that there is no place in the derivation where the actual solar flux is 
introduced. 
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Beard and Johnson in I960 discussed the interaction of a satellite 

8 

with the earth’s magnetic field. The induced potential gradient may 
be as high as 0.2 volts per meter which affects the distribution of the 
electron flow to the satellite surface and may also affect measurements 
of satellite pctential. Equilibrium potentials (at the midpoint of the 
satellite) are on the order of one volt negative for an electron temper- 
ature equivalent to 0.1 volts. Photo- and secondary emission were felt 
to be unimportant. 

Results of measuring the potential of Sputnik III have not been re- 
ported in detail. A summary of available statements indicates a neg- 
ative potential varying from -2 volts to -7 volts with altitude and with 
day -night conditions. 

Chopra in a review article found an expression for the potential of 
a body at rest: 


<£=- 




( 2 . 2 ) 


where m is the particle mass, the subscripts referring to ions or elec- 
trons.^ If photoemission predominates over ion collection, then 



e 



(2.3) 


where n e v e and n^ refer to the plasma electron flux and photoelectron 
flux respectively. At satellite velocities the ion flux is increased, but 
he believed that "at least in the outer parts of the terrestrial atmosphere 
and in the interplanetary space, the photoelectric effect is important." 

He observed that surface phenomena such as secondary emission are 
unimportant for particle impacts at ordinary gas temperatures. In a 


18 



paper in 1961, Beard and Johnson discussed ionospheric limitations on 
attainable satellite potentials. Higher negative than positive potentials 
are attainable by ejecting positive ions from a source in the satellite 
because of the limited mobility of environmental positive ions constitut- 
ing the return current. 

The probability distribution for charges on lunar dust grains was 

28 

considered by Grannis. However, he seems to have confused the 
cna.rge on a grain with the rate at which the grain acquires the charge 
and in addition did not take into account the effect of the grain's charge 
on its rate of charging. 

The latter criticism was also made of Grannis' paper by Walker 

who derived a different probability distribution for the charges on lunar 
29 

grains. In another paper Singer and Walker concluded that dust on 

the lunar surface has a charge proportional to the exposed surface 

area. ^ No large potential differences can exist because free electrons 

above the surface conduct currents efficiently. Dust ejected from the 

surface by meteors may become charged by the same processes that 

charge dust in interplanetary space. 

In a companion paper Singer and Walker calculated the screening 
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effect of photoelectrons on bodies in interplanetary space. Photo- 
emission current density was computed using Hinteregger 's results on 
the number of solar photons with energies greater than 8 eV. The yield 
was assumed to be 1 for want of better information. When applied to 
the lunar surface a potential of about plus 20 volts was obtained. 

Gdalevich has reported some results of electric field measrrements 
from rockets launched in 1957 and 1958.^’^ He found a field of 0.2 
to 3 volts /cm at the rocket surface corresponding for the most part to 
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a negative charge, although for portions of the trajectory the charge 
was positive. He also has derived expressions for the rocket potential 
similar to those of Gringauz and Zelikman and of Chopra. Imyanitov, 
in reporting results from more recent rockets, found rather high neg- 
ative potentials (several volts) which led him to assume that a consider- 

34 

able number of fast negative particles was present in the atmosphere. 

Sagalyn and others have also reported negative rocket potentials, from 

35 

-0.4 volts at 150 km to -1.7 volts at 450 km. 

Shen and Chopra have considered the problems accompanying 

36 

thermionic emission from a hot body in a plasma. Solutions were 
obtained for the potential resulting from the balance between thermionic 
emission and electron accretion. 

Rawer in a recent article has discussed the positive and negative 

37 

particle fluxes to a satellite. He has some good comments on the 
inductive effect of the earth's magnetic field, but states that the mag- 
netic field has no effect on the isotropy of the electron flux apart from 
the induced polarization. The effect of photoemission is discussed 
carefully with reliable values for the solar flux and yield. He points 
out that in regions of low electron density the satellite potential may 
be determined by strong emission lines such as Lyman alpha of H or 

He. + Finally, both Walker and Bettinger discuss equilibrium potential 
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in recent dissertations. * Walker finds a transcendental equation 


<*>=- 



ln 


/T^m. /fini 

e i i c 
l-0e/k T. 




(2.4) 


for the potential of a sphere at rest in a plasma. Bettinger computes 
the equilibrium potential for an insulated probe that emits electrons 
thermally. He neglects photoemission but assumes that there is a high 
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energy tail to the normal Maxwell-Boltzmann velocity distribution of 
electrons. The results of a rocket flight carrying the probe are pre- 
sented and discussed in terms of the number of high energy electrons 
necessary to obtain the observed values of -1 to -4.5 volts. 

2. Present Status of the Problem . It is apparent that a considerable 

amount of attention has been devoted in the literature to the problem of 

the equilibrium potential of a body in space. However, almost all of 

the treatments have been restricted to a consideration of two or perhaps 

three mechanisms that the author considers imports"' for his model; 

one exception to this has been the discussion of the charge on inter- 

2 

stellar grains, particularly Spitzer's work. Only a few experimental 
data have been obtained, and in no case have these measurements been 
analyzed in terms of the expected potential where all the possible 
charging mechanisms were evaluated for the specific vehicle that carried 
the experiment. 

In addition to Sputnik III, vehicle potentials have been measured on 

the satellites Explorer VIII, Ariel I and Explorer XVIi. Some prelim- 

40 

inary results have been reported from Explorer VIII. No data on 
satellite potential have been reported in the literature yet from Ariel I 
or Explorer XVII. 

The author believes that some of the cnarge acquisition mechanisms 

have not been treated sufficiently thoroughly. Only Rawer has attempted 

to evaluate the effects of photoemission by combining photoemission 

3 ^ 

yields as a function of wavelength with the solar spectrum. Our 
knowledge of the solar spectrum in the extreme ultra-violet has since 
been improved, and the author feels that fairly good quantitative esti- 
mates of the photoemiss.ion current can be made for certain materials. 


21 



The effects of energetic particle fluxes have been invoked as being 
of probable significance, but there have been no quantitative estimates 
based on known fluxes such as those of the Van Allen radiation belts. 

Such an estimate should be made, and secondary emission yields are 
known well enough for certain materials that this effect cculd also be 
included in the calculation. 

Another effect that has only been noted in passing and then dismissed 
is that of the magnetic field in restricting the direction of motion of 
environmental ions and electrons. It will be shown that this effect can 
be quite significant. 

Finally, there is no literature presently available should one want 
to make a quantitative estimate of the equilibrium potential for, say, a 
satellite or other body of certain dimensions and material under specific 
er. /ironmental conditions. It is hoped that this investigation will help 
to fill this gap. 
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CHAPTER III 


COLLECTION OF ELECTRONS AND IONS 

1 • General Considerations . Mechanisms of charge acquisition can be 
classified as charge collection or charge emission. The latter consists 
of processes such as photoemission, thermal emission, field emission 
and emission of alpha or beta particles from radioactive materials in 
the body. There may be a combination of collection and emission as 
when there is secondary emission of electrons upon incidence of ener- 
getic particles. It turns out that by far the most important processes 
are collection of environmental electrons and ions, which will be treated 
in this chapter, and photoemission and secondary emission, which will 
be discussed in Chapter IV. Other less important processes will be 
considered in Chapter V. 

At the outset a distinction may be made between the incidence of 
energetic particles and that of lower energy (thermal) particles in that 
only the latter are influenced by the charge on the target body. Hence, 
it is a straightforward calculation to determine the current to a body 
from energetic particles if the particle flux and directional distribution 
are known. The effective collection area for a unidirectional flux, for 
example, will simply be the cross-section of the body normal to the 
particle flux. The total effect of the energetic particles on the body's 
charge must of course take into account the amount of induced secondary 
emissions. This and the related problem of "sticking probabilities" 
for incident ions and electrons is considered in the next chapter. ' 
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Finally, a distinction should be made between large and small 
target bodies according to their dimensions compared to a Debye 
length, L : 



where e Q is the permittivity, k Boltzmann's constant, T the plasma 
temperature, n the electron or ion density and e the unit electron 
charge. The significance of the Debye length is that any shielding of a 
charged body by space charge in the surrounding plasma occurs in a 
distance on the order of a Debye length. Hence, for example, a spher- 
ical body that is small compared to the Debye length is effectively 
unshielded for many radii away from its center and the electric field 
is essentially coulomb. On the other hand, bodies that are large com- 
pared to L have their charge shielded in a small fraction of a radius 
away from the surface. Consequently, it is sometimes possible to 
treat the problem of particle attraction by assuming a neutral body with 
a slightly larger surface area. 

2. Electron Collection . The problem of electron and ion currents to 

a probe in a plasma has been the object of considerable attention in 

the literature, beginning with the work of Langmuir, Mott-Smith and 
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others three decades ago. ' In the upper atmosphere and in space 
the situation is somewhat simpler than in the laboratory in that there 
are no "wall effects" such as occur in laboratory vacuum systems, 
and the thermal plasma is probably more nearly Maxwellian. This is 
particularly true for electrons since their most probable thermal 
velocity is much larger than typical satellite or meteor velocities. 
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Velocities range from about 8 km/sec for a satellite near the earth to 
a maximum of 73 km/sec for a meterorite approaching the earth, 
whereas the electron thermal velocity is on the order of 200 km/sec 
for a temperature of 1500°K. Therefore it is realistic to assume that 
in a satellite centered co-ordinate system the electron velocity distri- 
bution is still Maxwellian. This is not true for ions which have thermal 
velocities on the order of 1 km/ sec — in general lower than typical 
satellite velocities. 

One other simplification is that at the altitudes which are considered 
here, collisions between particles are unimportant. The minimum 
perigee altitude for a satellite that is to have a lifetime of at least a 
few days is about 150 km. The mean free path at 150 km is about 50 
meters which is larger than most satellite dimensions and very large 
compared to a Debye length at that altitude. 

For a repulsive potential the electron current to a body will follow 
the Boltzmann relation. 


neaA <t>e/kT 

— F* e 

2 \/tt 


(3.2) 


where a is the most probable thermal velocity defined by 1/2 ma 2 = kT, 

and A is the surface area of the body. This equation is valid for any 

convex-shaped body of any size in contrast to the expressions for attrac- 

41 

tive potentials which depend strongly on the body's geometry. There 
are other phenomena, however, which may modify the effective collection 
area A , such as the effect of a magnetic field or the presence of space 
charge in the wake behind a satellite, as will be shown later. 
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For attractive potentials both the size and shape of the body are 
important. Spherical, cylindrical and planar geometries have been 
treated in the literature. It is convenient to restrict the discussion at 


this point to spherical bodies for several reasons: (1) It is clearly 
advantageous to approximate a complex- shaped body as a sphere for 
simplicity of treatment. (2) This is a good approximation for small iso- 
lated bodies where the far field will be a coulomb field regardless of 
the details of the body's shape. (3) Most satellites are roughly spher- 
ical in that the three axes are approximately equal. Exceptions come 
immediately to mind such as antennas or long booms, but they may be 
treated separately as special cases. 

A positively-charged sphere, then, whose radius is small compared 

to the Debye length will be surrounded by a coulomb field that is effec- 

40 

tively unshielded. The electron current to a body in such a case is 

I=n ^r( 1+ cf) (3 - 3 » 


For a body whose radius is comparable to the Debye length the 

variation of the potential through the sheath is important in determining 

the total collected current. This means that Poisson's equation must 

be solved in the sheath taking into account all the sources of space 

charge in addition to that due to the electrons alone. Ia general this is 

an extremely complicated problem requiring numerical procedures, but 

in certain simplified cases an approximate analytical expression may 

be obtained. Mott-Smith and Langmuir derived the following expression 

41 

for the current due to attracted particles in a spherical sheath: 
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(3.4) 


I - (4na 2 ) fl - (_ 8 - 2 " r2) e “r 2 e^/kT(8 2 “r 2 ) 

2/rr La 2 

where r is the radius of the body and a the radius of the sheath, assumed 
to be spherical and concentric with the body. This equation was derived 
by assuming a well-defined edge to the sheath so that the flux of elec- 
trons at the sheath surface is due to their random thermal motions. 

More recently it has been shown that the sheath boundary is not well- 

defined but that the electric field may penetrate "beyond" the sheath for 
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a considerable distance. However, Walker has found that the Mott- 
Smith- Langmuir equation may still be used if the following expressions 
are used to determine the sheath radius, 


a = r + t 


(3.5) 


and 


f = 0.83 

L 


(rHiT 


(3.6) 


where t is the "thickness" of the sheath around the body. Equation (3.4) 
reduces to equation (3.3) for the case when t >> r ; when t << r then 


I = (47ra 2 ) ^ (4 W r 2 ) (l +— ) (3.7) 

2At 2 /n \ >"/ 

i.e., the body may be taken to be neutral with a slightly larger ( 2 1 / r) 
surface area. 


3. Positive Ion Collection . The essential difference between the treat- 
ment of electron collection and ion collection is that in general the 
motion of the body through the plasma cannot be neglected. In the 
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special case where the body is at rest with respect to the plasma, the 
preceding equations derived for electrons are applicable to ion collec- 
tion with appropriate changes in the sign of the potential and in the 
quantities referring to particle characteristics. In succeeding equations 
the subscripts + or - will be used to denote reference to ions or elec- 
trons. These equations for the body at rest will also serve as checks 
on. the general ion current equations since the latter must reduce to the 
former for zero velocity. 

41 

Mott-Smith and Langmuir discussed the problem of a moving 

collector in a plasma but did not give equations for the current. Gringauz 
21 

and Zelikman showed that the ion current tc a moving sphere would 
decrease in a nearly linear manner with an increasing retarding poten- 
tial, approaching zero at a potential corresponding to the kinetic energy 
of the ions in the moving system, <f> = m + V 2 /2e, where V is the velocity 
of the sphere. The current is given approximately by 


I + = 7rr 2 n + eV 



(3.8) 


as long as the thermal motions of the ions can be disregarded. Two 
spherical ion traps, each consisting of an outer grid and an inner collec- 
tor biased to repel electrons, were flown on Sputnik III. Linear current- 

26 

voltage curves were obtained for repulsive potentials as predicted; 
effects of the ion thermal motions were discussed but a general equation 
was not given. 

Such an equation for repulsive potentials has been derived by 
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Hintereg^er and later by Kanal. 
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t nt \t 

I+ = i" n+eV 


a: 


2V 2 


U 2 

V 2 



+ 




(3.9) 


where 


U= [2e<^/m + ] 1/2 


(3.10) 


and 


erf x = 




(3.11) 


Interestingly, the current is independent of sheath size or variation of 
potential through the sheath as long as the electric field is radial. 
Gringauz showed that this is because the current is limited by angular 
momentum considerations, and is true as long as the effective radius 
of the collector, given by 


T efl = r 


1 2 e( M 1/2 

m + V l 


(3.12) 


is less than the sheath radius. This is always true for repulsive poten- 
tials since the bracket in (3.12) is less than unity. In addition, it will 
be true for attractive potentials until the effective radius exceeds the 
sheath radius, with the consequence that (3.8) may also be used for a 
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limited range of attractive potentials. Equation (3.9) reduces to (3.8) 
as the ratio a + /V approaches zero, and it reduces to Equation (3.2) as 
a + /V approaches infinity. 

An analytic expression for the general case of the ion current to a 
moving body in an attractive field has not yet been obtained. The prob- 
lem is extremely complex because it involves a simultaneous solution 
of Poisson's equation and calculation of the ion trajectories. Both 
Kanal^ and Walker^ have made certain simplifying assumptions in 
order to obtain useful solutions. Kanal used a model for the sheath in 

I 

which ions enter the sheath with zero initial velocity and obtained two 
equations for the ion current by means of which the unknown sheath 
radius was eliminated graphically. Walker integrated the ion trajectories 
numerically, assuming a spherical sheath edge. The shape of the collec- 
tor then depends on the results of the computations, and is not in gen- 
eral spherical. For the case where the collector is nearly spherical, 
a comparison of Kanal's and Walker's currents shows a disagreement 
by more than an order of magnitude. 

The assumption of a spherical sheath in front of the moving body is 

probably quite realistic as has been shown by Al'pert, Gurevic and 
.. 46 

Pitaevskij. If the sheath edge were sharp it would be possible to use 
Equation (3.9) for the range of potentials where r eff < a , the sheath 
radius. The same equation with <p set equal to zero, and with r re- 
placed by a could be used for the case when r eff > a . Unfortunately, 
the sheath edge is not sharp as Walker and others have pointed out. 
However, the precise position of the edge of the sheath will not matter 

greatly if the sheath thickness t = a - r «r. In this case the current 

35 

will be the same as that derived by Sagalyn et al. 
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(3.13) 


I + = 7ra 2 n + eV 


a* 


a 


2V 2 j 


erf 


+ *(V/a + )' 


a. 


V A 


In the case when the body is so small compared to the Debye length 

that the field may be assumed to be a coulomb field, the current for an 

45 

attractive potential has been derived by Kanal. He finds 



The current in this case reduces to that in (3.3) as the velocity of the 
body approaches zero, and it reduces to (3.13) for zero potential. 

Because of the importance of the ion current to a negatively charged 
moving sphere for satellite potential calculations, a program to com- 
pute this current has been developed. Poisson's equation has been 
solved with the assumption of spherical symmetry and with the space 
charge given by 


p = n + e (1 - e ^ e/kT ) (3.15) 

i.e., the ion density is constant and the electron density is described 
by the Boltzmann factor, the same assumptions used by Jastrow and 

7 

Pearse. The ion current may then be computed, without integrating 

the ion trajectories, by using Walker's classification of trajectories as 

38 

either periastron or pericritical. Details of the analysis are given 
in the appendix. Some typical results are given in Figure 4 in the form 
of current versus satellite potential curves. The linear relationship at 
lower voltages with a saturation effect as the effective radius of Equa- 
tion (3.12) exceeds the sheath radius is apparent. The thermal motions 
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EQUATION (3.8) 
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Figure 4. Ion Current to a Large, Negative, Moving Body 




of the ions have been neglected, but the temperature of the plasma is 
taken into account through the Debye length. 

The accompanying table summarizes the different situations that 
have been discussed and shows the appropriate equation to use in each 
case. 

Table I 

Size of Body: Small Intermediate Large 

Body at Rest, V/a < < 1 


Attractive 

Field 

Eq. (3.3) 

Eq. ( 3.4) with 
(3.5) and (3.6) 

Eq. (3.7) 

Repulsive 

Field 

Eq. (3.2) 

Eq. (3.2) 

Eq. (3.2) 



Body Moving, 

V/a/0 


Attractive 

Field 

Eq. (3.14) 

Fig. 4 

Eq. (3.13) 
or Fig. 4 

Repulsive 

Field 

Eq. (3.9) 

Eq. (3.9) 

Eq. (3.9) 


4. Effects of the Satellite Wake : A body moving rapidly through the 
atmosphere (in comparison with the ion thermal velocity) will have a 
rarified region behind it in the shape of a cone, as shown in Figure 5, 
with a half-angle given approximately by tan 0 - a + /V. The shape of 
the wake and the ion and electron and potential distribution in the wake 

4f 

have been discussed in great detail by Al'pert, Gurevic and Pitaevskij. 
The potential distribution in the wake does not appreciably affect the 
total ion current to the body because most of the ion current is incident 
on the front half. However, this is not the case for the electron current. 
Electrons will diffuse into the wake where because of the absence of 
ions there will be a net negative space charge. Al'pert et al. show that 
the potential in the wake will in general be more negative than that of 
the body, reaching an extremum of approximately - 2kT/e In (r/L) as 
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shown by the dotted line. Consequently, electrons with energies less 
than this will be unable to cross this potential barrier to reach the rear 
half of the body, and in place of Equation (3.2) one should use 


I_ = 


n_eu_ (27/t 2 ) 

2 V TT 



(3.16) 


as long as 4 > s , the satellite potential, is more positive than - 2kT/e ln(r/L), 
but less than zero. If the satellite potential is positive, then Equation 
(3.7) should be used for the front half of the satellite, and 


I 


n_ ea_ (2ttt 2 ) 

2/77 



( 3 . 17 ) 


for the rear half of the satellite. 
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CHAPTER IV 


PHOTOEMISSION AND SECONDARY EMISSION 

1* Photoemis sion . Two steps are involved in estimating the effects of 
both photoemission and secondary emission on the equilibrium potential 
of a body. Only the first step is necessary if the body is negative, 
namely the calculation of the total emission current which depends on 
the energy spectrum of the primaries and the secondary electron yield 
as a function of primary energy. Since the body is negative all the sec- 
ondaries may be assumed to escape, and it is unnecessary to know the 
energy spectrum of the secondary electrons. 

If the body is positive not all of the secondary electrons will escape, 
and it is in this case necessary to know the energy spectrum of the sec- 
ondaries. If the body is small (compared to a Debye length) then all 
electrons emitted with energies greater than that corresponding to the 
potential difference between the body and the environment will escape. 

If the body is large the equipotential surfaces will be approximately 
planar, and the condition for escape is that the directed kinetic energy 
of the electrons normal to the surface must be greater than4> s e, where 
$ s is the satellite potential. In general, for a sphere we have 


E> I V(r) + 


= [V'(r>], 


( 4 : 1 ) 


as the condition for escape, where E is the total energy (kinetic plus 
potential) upon emission, V(r) is the potential outside the sphere, and 
p is the angular momentum of the emitted electron. 
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The importance of photoemission in providing a charging mechanism 
for bodies in the upper atmosphere and in interplanetary space has 
emei 0 ed as a result of rather recent experimental data in two areas of 
investigation. First has been the demonstration of rather large photo- 
electric yields for metals in the extreme ultra-violet range of wave- 
lengths— ranging up to peak values greater than 10% in some cases 
compared with numbers like 10“ 4 to 10" s for yields in the near UV 
(> 2000.X.) . And second has been the discovery of considerable energy 
in this same region of the solar spectrum. Figure 6 shows the solar 
spectrum from 100 to lO.OOoJL The ordinate is the number of photons 

per cm 2 per second per Angstrom interval. The data at wavelengths 

o 47 

below 1775A was obtained by Hinteregger, Hall and Schmidtke; that 

O 4g o 

above 1775A is from Nawrocki and Papa. The Lyman-a line at 1216A 
is particularly important in its contribution to photoemission. 

To demonstrate the effects of both photoemission and secondary 
emission upon the problem of equilibrium potential, two metals have 
been chosen: tungsten with a work function between 4 and 5 electron 
volts and aluminum with a lower work function between 3 and 4 volts. 

The choice of these materials is partly dictated by the fact that data is 
available for both secondary and photoemission from both. Also, alumi- 
num is used extensively in spacecraft construction, and tungsten has 
been used frequently in experiment sensors so that a comparison of 
laboratory and flight results is possible. 

Photoelectric yields for tungsten in the ultra-violet wavelengths 

have been measured both in the laboratory and on rocket flights in the 

upper atmosphere. The yields as a function of wavelength are shown 

in Figure 6 for both clean and dirty tungsten: "dirty" means an untreated 

surface, and "clean" means the surface was heated at a temperature 
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Figure 6. The Solar Spectrum and Photoemission Yields 
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greater than 1000°C in a vacuum of 10” 5 torr until yield reproducibility 
was established. Hinteregger 's data was obtained for an untreated sur- 
face; he states that the yields are reproducible even after exposure to 
air. 

The solid line drawn through Hinteregger and Watanabe's experi- 
mental points to higher wavelengths is a theoretical curve for the yield 
according to the Fowler -DuBridge theory, with a long wavelength cut-off 

O 

at 1900A, corresponding to a work function of 6.5 eV determined by 

56 57 

Fowler's method. ’ The theoretical curve through Rentschler's 

O 

points corresponds to a long wavelength cut-off of 2690A, or 4.6 eV, in 

58 

good agreement with Warner's results. The difference in the curves 
is probably due to the state of the metal surface, although it is not clear 
that the Fowler -DuBridge theory can be applied at wavelengths as low 
as 1700A. 

When the product of the yield and solar flux is integrated over the 
spectrum, a total current of 2.1 X 1CT 9 amp/cm 2 and 8.1 X 10“ 9 amp/ 
cm 2 is obtained, depending on which of the two curves is used at the 
longer wavelengths. These values bracket experimental photoemission 
current densities Lora tungsten of 3.9 X 10“ 9 amp/cm 2 obtained by 

Hinteregger^ and 5 X 1(T 9 amp/cm 2 obtained on Explorer VIII.^® 

0 

The yields for aluminum below 2500A are laboratory measurements 
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by Suhrman and Pietrzyk fitted to the following law by Rawer: 


r *.2 x KT 3 (<D X - $> w ) 3 , 4> x < 8 volts 

Y = 0.077, > 8 volts 


(4.2) - 


where is the work function of 4eV and <t>^ is the energy of the cor- 
responding wavelength. These values join nicely with experimental 
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data by de Laszlo^ above 2500A. These large yields, which are a con- 
sequence of aluminum's lower work function, result in a computed pho- 
toemission current density of 2.5 X 10" 7 amp/cm 2 above the atmosphere- 
two orders of magnitude greater than that for tungsten. Such a large 
photocurrent would certainly be important for satellite potentials because 
of this metal's wide use in spacecraft construction. Consequently, an 
experiment was designed to measure the photoemissior- current from 
aluminum and was flown on a Nike-Apache rocF to an altitude of 193 
km on Dec. 16, 1964 at White Sands, IV. M. The maximum measured 
current at an altitude of 160 km was 2.3 X 10“ 9 amp/cm 2 . An extrapola- 
tion to the top of the atmosphere yields a value of 3 X 10~ 9 amp/cm 2 , in 
fair agreement with the result for tungsten rather than the computed 
value for aluminum. Probably the discrepancy is due to the state of the 
aluminum surface causing an increased work function — de Laszlo's 
results were for a carefully prepared and outgassed surface. The ex- 
perimental photocurrent will be used in the calculations of satellite 
potential. 

When a satellite or other body is positive, the energy distribution 

of the photoelectrons is needed. An experimental cur rent -voltage curve 

far photoelectrons emitted from a positive plane tungsten surface has 

50 

been obtained by Hinteregger, Damon and Hall, and is reproduced in 
Figure 7. For small bodies where the angular momentum distribution 
of the electrons is needed to compute the number that escape, it will be 
assumed that they are emitted with a cosine distribution. 

2. Secondary Emission of Electrons Upon Electron Impact. The phe- 
nomenon of secondary electron emission upon electron impact has been 
studied extensively and is fairly well understood according to recent 
reviews of the subject.^’^’^* »k2>63 The electrons emitted from a 
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target surface include, however, reflected and back-diffused primaries 


as well as true secondary electrons. The total yield, 8 , defined as I t /I p , 
the ratio of emitted target current to incident primary current, may be 
written as 


8 = S p (1 _ r) + ,■ + 77(1 + 8 p /3) (4.3) 

Here 8 p is the "true" yield of secondaries due to the primaries, r is the 
reflection coefficient (1-r is the "sticking factor" for electrons), 17 the 
back-diffusion coefficient and fi an efficiency factor describing the 
increased efficiency with which the back-diffused electrons produce 
their own secondaries. 

The shape of the true yield curve as a function of primary energy 
below a few keV is in effect a universal curve for metals when normal- 
ized to the maximum yield and the primary energy at the maximum. 
Figure 8 shows the yields for A1 and W as a function of primary energy. 
What is plotted in Figure 8 is the effective yield given by 8 p (l + 7 ?/^)> 
since measurements of the yield generally have been corrected for the 
back-diffused primaries only, with no distinction being made between 
slow secondary electrons caused by the incident as against the back- 

scattered primaries. The yields have been extrapolated to zero at 5 eV 
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according to the observations of Shulmanand Myakinin, andHarrower 

that true secondaries occur only above primary energies of 4-7 eV. The 

maximum yield occurs approximately at an energy where the primaries 

are stopped at a depth in the metal corresponding to the depth from 

which secondaries can escape. At higher primary energies the yield 

decreases in accordance with the rate atwhich the primaries lose energy 

in the escape zone; and it has been shown by Kanter^ for energies 
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between 1 and 20 keV and by Schultz and Pomerantz ' at energies up 
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to 1.6 MeV that the energy loss rate is in agret-.ient with Beihe's 
stopping -power formula: 


5 p r k S = ¥ ln <E/Cj) (4 ‘ 4) 

The reflection coefficient, r , is significant only at low primary 
energies below about 10 eV. In general, r is of the order of 0.05 at zero 
primary energy and decreases with increasing energy according tr> the 
relation 


r = (4.5) 

16 (E + W ) 3 +W 3 

where W is the sum of the Fermi energy and work function of the metal, 

Z Q 

and E the primary energy, both measured in Rydberg units ( 13.54 eV). 
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Guth and Mullin have found r to be 0.05 for tungsten n^ar zero volts. 
The use of Equation (4.5) with a work function of 3.5 eV and a Fermi 
energy of 5.6 eV yields a reflection coefficient of 0.04 for A1 at zero 
volts . 

Figure 9 shows the back-diffusion coefficient 17 for A1 and W as a 
function of primary energy. The curves have been extrapolated to 
zero at 100 eV because the average energy, E/E of back-diffused elec- 
trons is 0.50 for A1 and 0.60 for W in terms of the primary energy, 
whereas 50 eV is usually taken to be the energy distinguishing back- 
diffused electrons from true secondaries. 

Figure 10 shows the total yield 8 for Ai and W as a function of 
primary energy, obtained by combining the appropriate values for the 
various coefficients in accordance with Equation (4.3). 
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Figure 10. Total Secondary Electron Yield for Tungsten and Aluminum 



It should be added that above 8 p max the yield is dependent upon the 
angle of incidence of the primaries. At energies near the maximum 
yield it is given by 


In (8/8 0 > =C (1 - cos 9 ) (4.6) 
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where C depends on the energy of the primaries. At higher energies 
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the yield is proportional to sec 9 . An effective yield for an isotropic 
flux of primaries may be defined by 


8 


elf 


J .ir/2 

0 


8(6f) sin 9 cos 9 d 9 


(4.7) 


This reduces to 8 = 8. for the 8 = 8. sec 9 law. 

e f f 0 O 

All the secondary electrons emitted from a spherical body in space 
will escape if the satellite potential is negative. For positive potentials 
it is necessary to know the energy distribution of the secondaries; how- 
ever, for reasonable potentials— i.e. less than +50 volts— only the energy 
distribution of the true secondaries is needed since the back-diffused 
primaries will still escape. 

Schultz and Pomerantz state that "the energy spectra of secondary 
electrons emitted from metals bombarded by relativistic electrons are 
practically identical with those measured at very much lower primary 

/ *7 

energies by Kollath." Figure 11 gives this differential energy spec- 
trum. There is abundant evidence that the. angular distribution of the 
secondaries follows a cosine law, hence it is a straightforward pro- 
cedure to compute the number that escape from a body with the use of 
Equation (4.1). 

3. Secondary Emission of Electrons Upon Ion Impact . The predominant 
positive ions in the upper ionosphere are 0+ from 150 km to approximately 
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igure 11. Energy Spectr 





10 3 km and H + at higher altitudes with a more or less thin belt of He + 
between. These are thermal ions, with energies on the order of a frac- 
tion of an eV. In addition, there are belts of energetic protons, with 
energies up into the 100 MeV range in the radiation belts. 

Below about one keV the yield for secondary electron emission 
from ion impact is very nearly independent of the ion's kinetic energy. 
This is because the ion is neutralized directly to the ground state as it 
approaches very near to the metal surface. The number of excited 
electrons depends on the available potential energy after neutralization 
which is determined by the ionization potential, <£. , of the incident atom 
and the work function of the target metal , as illustrated in Figure 12. 

When a conduction electron is captured by the incident ion, it makes 
available a maximum energy of <£. - 0 W • At least 4> m of this must be 
used to free another electron from the metal so that the condition for 
secondary emission by this mechanism is that </>. >2 <£ w . It is apparent 
from the figure that the yields for various ions incident on several 
metals depend primarily on the difference 0. _ 2 This relationship 
is used to estimate the yields in the following table with the exception 
of He + on W which has been measured by Hagstrum. 


Table II 



Al, = 3.5 eV 

W, <f> = 4.6 eV 

H + , 4>. = 13.5 eV 

0.086 

0.041 

He+, 0. = 24.6 eV 

0.38 

0.295 

0+,^ = 13.6 eV 

0.088 

0.045 


It should be added that the yield is quite dependent on the condition 
of the metal surface. Hagstrum's data is for atomically clean surfaces, 
whereas the platinum in Parker's measurement may have had some 
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Figure 12. Secondary Elect 
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residual impurity even though the measurement was made only 25 sec- 
ondr after flashing the surface. The He* + ion probably liberates elec- 
trons through the two-step resonance capture process rather than by 
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the Auger process. 

Figure 13 shows the yields for H* on A1 and W from 1 keV up to 

about 10 MeV. The curves have been extrapolated to the value of Table 

II at 1 keV. They show a broad maximum at about 100 keV converging 

to a single curve independent of the target material at higher energies. 

The variation of the ion yield as a function of angle of incidence of 
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the primaries has been measured by Oliphant and by Allen. Their 

results show that the yield is proportional to sec 9 . Thus the effective 

yield for an isotropic flux of monoenergetic ions incident on a surface 

will be unchanged if the flux is given in units of no/cm 2 -sec-steradian. 

Reflection coefficients for He + incident on clean and contaminated 

tungsten have been measured by Hagstrum to be 0.0017 and 0.00043 
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respect /elyand were found to be fairly insensitive to the ion energy. 
Although reflection coefficients for H + and 0 + on tungsten or aluminum 
are not available, they are probably somewhat larger but still less than 
0.02. Hagstrum suggests that these values are representative of the 
reflection of ions whose ionization energy is large compared to twice 
the work function of the solid. It is concluded that ion reflection is 
unimportant as far as charging effects on a body in space are concerned. 

The energy distributions of secondary electrons are quite dependent 
on the state of the surface but in general peak near 2 eV with a Max- 
wellian shape. For practical reasons Figure 11 will be used when 
needed to compute the escaping electron current for positive satellite 
potentials. It will be assumed that the angular distribution of the 
secondary electrons is a cosine function. 
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CHAPTER V 


OTHER CHARGING MECHANISMS 

1. Discharge Tim e for a Body in a Plasma. It was stated in Chapter 
III that the most important mechanisms for charging a body are the 
collection of environmental electrons and ions and the secondary proc- 
esses treated in Chapter IV. It is necessary to consider other possible 
mechanisms to see under what conditions this assumption is justified. 

The concept of the discharge time r for a body in a plasma is use- 
ful in assessing the relative importance of a proposed charging mecha- 
nism in comparison with ion and electron collection. Suppose that a 
body at an equilibrium potential determined by a b iance between ion 
and electron accretion from the environmental plasma suddenly acquires 
an additional amount of charge by some other mechanism which acts 
only momentarily. The body will return to its previous potential in a 
time ou the order of r seconds. If the frequency of occurrence of this 
other mechanism is small compar ed with the quantity r" 1 , then its 
effect on the average equilibrium potential may be ignored. 

As an example, consider a small body at a negative potential where 
the ion and electron currents are given by Equations (3.2) and (3.3) 
respectively. The chan ge in the charging current to the body due to a 
sudden small change in potential is on the order of 

. _ n e a A e _ , . 

A I ^ — — A cp (5*1) 

2 /tt kT 
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Hence, we may write 


dQ “ n e a + A e Q 0 

dt " 2 ]/tt CkT t 


where Q is the excess charge over that at equilibrium and C is the body’s 
capacitance. Thus we find that the discharge time t for a small sphere 
of radius r is given by 


2 /7ir CkT 

T r 

n e 2 a + A 



(5.3) 


where r~ ! is the ion plasma frequency, 1/27? /n e 2 A 0 »» + , and t 0 is a 

characteristic time describing how long it takes an ion to travel a dis- 
tance comparable to the dimensions ~f the body. Representative values 
for r for a 10 micron radius body are 2 X 10 -3 sec in the ionosphere 
and 2 X 10 + 3 sec in interplanetary space. Discharge times for larger 
bodies will be proportionately smaller although sheath effects will then 

complicate the definition of r. 

1 7 

2. Cosmic Rays. Spitzer has shown that cosmic rays have a negli- 
gible influence on the rate of charging of small dust grains in inter- 
stellar space. The flux of primary cosmic rays outside the earth's 
atmosphere is about 1 cm~ 2 sec -1 , whereas the flux of ions from the 
plasma in interplanetary space is at least 10 s cm -2 sec -1 . The ioni- 
zation produced in a particle along the path of the cosmic ray will be 
on the order of or less than 10 6 electrons/cm, and only those produced 
within a depth of 10 -7 cm from the surface will be able to escape. Hence, 
a cosmic ray particle is no more effective than a single low energy 
environmental electron in charging a small particle. 
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In large bodies the effects of showers in increasing the charging 
efficiency must be considered. If the body dimensions are large com- 
pared to the interaction length for star production, then both the pri- 
mary and the products of the stars will come to rest in the body. One 
would expect a maximum in the charging efficiency at a body dimension 
on the order of a few interaction lengths. There is an analogy between 
this picture and the situation in the earth's atmosphere where the flux 
of cosmic ray particles reaches a maximum at a depth of about 60 gm/ 

cm 2 , or one interaction length. The corresponding distance is on the 

♦ 

order of 0.1 to 1 meter for appropriate densities. However, the flux 
in the atmosphere at the Pfotzer maximum is only double that outside 
the atmosphere; hence it may safely be concluded that the increase in 
charging efficiency is less than an order of magnitude, and the charg- 
ing effect of cosmic rays may be disregarded. 

3. Radioactivity . Radioactive material in a body in space constitutes 
a charging mechanism both through the escape of emitted charged pri- 
maries from the radioactive nuclei and also through the escape of sec- 
ondary electrons excited by the primary in its passage to the surface. 

To compete with the minimum expected ion flux in interplanetary 
space of 10 s cm" 2 sec" 1 used in the preceding section, a surface ac- 
tivity of 3 p -curies/cm 2 is required if only the primaries are taken 
into account. The efficiency with which secondaries can be produced ' 
will be of the same order of magnitude as the yields for secondary 
electrons discussed in Chapter IV. Thus /? radiation has a low effi- 
ciency near 1% until it has slowed down to less than one KV, where the 
efficiency peaks at values somewhat greater than unity. Proton and 
a -emission efficiencies will have a maximum value near ten at ener- 
gies near 100 KV. 
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This line of reasoning is supported by looking at the specific ioni- 
zation rates for these radiations. Peak ionization rates are less than 
10 ion-pairs /cm for both protons and a-particles in tungsten and occur 
near 100 KV, similar to the secondary yields in Figure 13. These spe- 
cific ionization rates were computed from stopping cross-sections per 
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atom, using 30 eV as the energy required to form an ion pair. The 
peak in production for electrons occurs at a few hundred eV. Ionization 
cross-sections for electrons in gases show a maximum less than about 
10- 15 cni 2 at these energies, which, if applied to a solid material such 
as tungsten, yield a specific ionization somewhat under 10 8 ion-pairs/ 
cm 60 . Since the secondary electrons will have energies below about 

O 

20 eV, only those excited within about 20A of the surface will escape. 
This yields an upper limit of 20 for the efficiency of a primary in pro- 
ducing secondaries which can escape from the body. 

The amount of radioactive material in bodies in space has been 
studied extensively from the viewpoint of determining the ages of mete- 
orites. Typical values for the induced activity as a result of cosmic 

ray bombardment are on the order of a few hundred disintegrations 
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per minute (dpm) per kg of material, or 10 -n curies /kg. Similar 
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activities have been measured on sections of recovered satellites. 

It is obvious that the charging effects of such low quantities of radio- 
active material are entirely negligible. 

Satellites sometimes carry quantities of radioactive material in 
conjunction with certain types of experiments, or as a power source; 
Such sources are usually well shielded but should still be considered 
as potential charging mechanisms. Clearly, each such source must be 
evaluated individually.' 
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4. Thermionic Emission . Shen and Chopra^ have discussed the effect 
upon the equilibrium potential of thermionic emission of electrons from 
bodies in space. They found that the temperature of a metallic body 
with a work function of 3.8 eV must exceed 700°K before it could have 
any significant effect. They assumed that the emission current was 
balanced by electron collection from the surrounding plasma charac- 
terized by a density of 10 3 /cm 3 and a temperature of about 1000°K. 

If a work function of 3.5 eV is assumed, a surface temperature of 
about 800°K is required for an electron emission of 10 s cm -2 sec -1 to 
compete with the minimum expected ion flux. Such a high temperature 
is likely only in special situations, such as when a meteor or space 
vehicle enters the earth’s atmosphere below 100 km at high speeds, or 
when a body approaches sufficiently close to the sun. A body with a 
typical albedo for meteors of about 0.4, whose temperature is deter- 
mined by a balance between solar flux and black-body radiation, will 
have a temperature of 800°K only when it approaches a distance of 0.19 
astronomical units— half the distance of Mercury’s orbit from the sun. 
The effects of thermionic emission will not be pursued further in this 
investigation. 

5. Field Emission . Spitzer and Savedof;:' have pointed out that for 
small dust grains field emission o£ electrons will limit the potential 
when it is negative. The onset of field emission occurs at surface 
field values between 10 6 and 10 7 volts/cm, the lower field correspond- 
ing to an emission flux of 10 s cm -2 sec~ l for a work function of 3.5 eV. 
Since the surface field of a small spherical particle is</>/r, the potent 

! *1 is limited to negative values below approximately 10 6 r. Thus, field 
emission is important only for dust grains of radius ~ 10“ 5 cm (0.1 
micron) and less. 



6. Collisions with Dust Grains . One square cm of satellite surface will 

be hit by a dust grain with a radius larger than 0.3 microns approximately 
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once every 50 seconds in the vicinity of the earth. This impact fre- 
quency decreases rapidly for larger grains; it probably also decreases 
by at least three orders of magnitude as one proceeds from the vicinity 
of the earth into interplanetary space. The value of 0.3 micron was 
chosen because smaller particles cannot remain at the earth's distance 
from the sun on account of the solar radiation pressure. 

In spite of the fact that one impact will produce many free electrons, 
the time between impacts is so large compared with the discharge time 
due to the plasma that the charging effects of such encounters may be 
completely disregarded. 

7. The Effect of Radio Frequency Fields . Early speculations**’^® that 

the rf fields around telemetry and other antennas on satellites might 

influence the satellite potential by a rectifying effect on currents from 

the plasma have been put on a quantitative basis by a group of Japanese 
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investigators working on radio frequency probes. ’ (See also recent 
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articles by Crawford and Harp ' on the theory of the rf resonance 
probe.) 

At rf frequencies below the plasma frequency the electron current 
density to a negative body is given by 

i = ioM v/v -) ( 5 * 4 > 

Here j 0 is the current density to the probe with no rf field, I 0 is the 
modified Bessel function of the first kind, v the amplitude of the rf volt- 
age and V_ the equivalent electron temperature. The current density is 
independent of frequency until close to the plasma frequency where it rises 

to a maximum and then falls to the value j fl at higher frequencies. 

60 



The Bessel function in (5.4) was obtained by integrating Equation 
(3.2) over one period of the rf voltage. Thus one would expect (5.4) to 
hold only at frequencies low enough that the electron's transit time from 
outside the sheath to the probe is small compared to the rf period. At 
positive potentials the effect would also be computed for low frequencies 
by averaging the appropriate equation from Table 1 over a period. 

The maximum in the current occurs close to but generally at a 
lower frequency than the plasma frequency. The resonant frequency 
depends primarily on the geometry of tht ^robe. The height of the cur- 
rent maximum depends both upon the geometry and upon the effects of 
damping by the plasma. Coliisional damping predominates if the neutral 

density is high enough; otherwise there is an "rf phase -mixing technique 
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akin to Landau damping." The former mechanism has been studied 
to some extent, but very little is known about the latter, which would be 
the prevailing mechanism in the upper atmosphere. The effect will be 
discussed further when the Explorer VIII data is presented. 
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CHAPTER VI 


MAGNETIC FIELD EFFECTS 

—4 _+ Q 

1* The Induced VxB Potential Gradient . Beard and Johnson pointed 

out that a satellite moving across the magnetic field of the earth would 

experience an induced potential gradient of as much as 0.2 volts/meter. 

This effect was observed on Explorer VIII through its effect on the elec- 

40 

tron current to a plasma probe. The phenomenon may be described 
as follows: a conducting body moving across the magnetic field is 
polarized so that in a coordinate system moving with the body the elec- 
tric field due to the induced polarization charge exactly cancels the 
induction field in the interior of the body. The polarization charge, in 
turn, is the source of a real field external to and carried along with the 
body which depends on the geometry of the body (and the environmental 
plasma), but whose effect on the plasma (i.e. on particle collection or 

on the sheath) may be described by saying that the potential of the sur- 

^ — > 

face of the body varies linearly with distance in the VxB direction. 

The best way to see this is to make a Lorentz transformation of 
the magnetic field to the moving satellite reference frame. A uniform 
electric field is obtained given by V * B in MKS units. The effect of 
this uniform electric field on a conducting body can then be calculated 
by the usual methods of electrostatics. For example, a sphere in a 
uniform field gives rise to a potential distribution described by a dipole 
term plus the uniform '‘.eld term. The effective external electric field 
is then obtained by subtracting the uniform field since the effect of the 
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magnetic field on the charged particles in the plasma can best be de- 
scribed by using the Lorentz force. Hence, the effective potential on 

— • *4 

the conducting surface is a linearly increasing function in the VxB 
direction. 

This effect is important primarily for the collection of electrons 
which depends exponentially on the satellite potential if it is negative. 
Beard and Johnson computed the effect on electron collection to a rec- 

g 

tangular parallelepiped satellite. The electron current to a negative 
sphere moving in a magnetic field may be computed as follows: in a 
spherical coordinate system centered in the sphere the surface potential 
is given by 

<£ = 0 o +|VxB|Rcos0<O (6.1) 

•—4 — + 

where the z -axis is taken in the VxB direction. An element of area 
given by 277 R 2 sin 6dd is at the potential given by (6.1), so that the cur- 
rent to the satellite is, from Equation (3.2), 


1. = -^ P e < *o* |,X3|RC "' ,) <e/kT) 2w R=sin«de 

2 Stt J 0 


(6*2) 


We obtain 


I 0 sinh (e |VxB| R/kT) 

(e | V x B | R/kT) 

where I 0 is the electron current that would be collected by the sphere 
if it were at the uniform potential 4> 0 . 

The practical consequence of this induction effect is that for large 
spacecraft one end of the structure is "pinned" close to zero volts or 
slightly positive, with the other end becoming relatively more negative. 



This is because the positive current to the body is limited, whereas the 
negative electron current increases more strongly. Thus, in Equation 
(6.3), as the V * B term increases, the potential 0 O contained in I 0 must 
become more negative to maintain 1^ relatively constant and equal to 

*4 * 

2. Effect of a Magnetic Field on the Direct Collection of Particles. In 
addition to the induction effect caused by the motion of a body, a mag- 
netic field can also affect the collection of charged particles by restrict 
ing. their direction of incidence. Consider, for example, a cylinder with 
its axis parallel to the magnetic field. Electrons approach the cylinder 
by spiralling along the magnetic field with a certain radius of gyration, 
p. For an electron to be collected at a given point on the side of the 
cylinder at least two conditions must be met: first, the distance of the 
magnetic field line about which the electron moves from the cylinder 
axis must be less than R + p, but greater than the absolute value |R _ p\, 
where R is the radius of the cylinder; and second, the distance of the 
point from the end of the cylinder must be less than the "wavelength" 
of the spiral. It is apparent that the current to the side of such a cyl- 
inder must decrease as the distance from either end increases. 

This effect is important whenever p is small compared to the size 
of a body, and consequently is important for electron collection by 
rockets and satellites, since p is on the order of cm in the ionosphere. 
Because of the fact that this effect on electron collection in the upper 
atmosphere has not been treated quantitatively in the literature, a cal- 
culation is made here for both cylindrical and spherical bodies with the 
assumption of zero potential difference between the body and the plasma 
The results for a cylinder will be applicable to probes mounted on the 
side of a rocket; indeed, there is evidence for anomalously low electron 
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currents to such probes on rockets flown at Fort Churchill, Canada, 
where the rocket axis was within a few degrees of the magnetic field 
direction. 

A Maxwell -Boltzmann velocity distribution for electrons may be 
written as follows, using cylindrical coordinates in velocity space: 


f (v , v T ) = n. 


/ m_ \3/2 _ 

\ 277 kT/ 6 


2kT 


< V ’ + V J> 


(6.4) 


where v z is the component of velocity along the field line, and v T is the 
transverse component. The radius of gyration is given by 

P- m_v T /Be (6.5) 


and the wavelength by 

X = 277m_v 2 /Be (6.6) 

At the plane z = 0 defining the front face of the cylinder, the current 
density of incident electrons is given by 


J z = JJ e v z v t^ v t ^ v t (6*7) 

where 2v v T dv T dv r is the volume element in velocity space. The geom- 
etry of the problem is described in Figure 14, where r is the distance' 
of the guiding center field-line from the cylinder axis, and 0 is the 
phase angle of the electron in its circular motion defined with respect 
to the 0=0 line as shown in the figure. All electrons with r < R-£>,^ 
where p < R , will be collected on the front face independent of its phase 
angle 0 at the plane z = 0. To be collected on the side of the cylinder 
the condition R~p<r<R 4 p must be met, and the phase angle must be 
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between 0 and 4> c as it reaches the plane z = 0. The particle impacts 
the cylinder when <p reaches 4> c at a distance z along the cylinder that 
is proportional to the phase « gle (<t> c -<£) and the velocity v s : 


z 


angular velocity 2 n c 


( 6 . 8 ) 


where 


<f> c = 2w-2 arc cos j (6.9) 

Obviously, if this value of z exceeds the cylinder length L, the particle 
will miss the cylinder completely. It follows that if the above condi- 
tions are satisfied for a given r, p and X, that all the particles with an 
initial phase angle <f> and with wavelengths between X and X v dX will be 
collected in the surface element 27rRdz, where 


dz = ~ (<fi c - 0) dX 


( 6 . 10 ) 


The current per unit length to the cylinder will then be given by 


dl _ J 2tttAt ( d<ft \ 

dz J dz \2nJ 


or. 


dl 

dz 


2n. ea_z .. r^c 

Jj 





dpdr d <f) 


( 6 . 11 ) 
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where a. has been defined previously in Chapter III, and where p* = 
m_a_/Be. The area ot integration for the variables p and r is the shaded 
area in Figure 14, bounded oy the lines • R + r, p R - r and p - r - R. 
This may be integrated with res pert to .* to obtain 


dl 

dz 


n_ ea_ 




dp dr 


( 6 . 12 ) 


So far we have been considering electrons incident at one end of the 
cylinder only. Electrons arriving from the opposite end will constitute 
a current of the same form but with z replaced by (L -z ). Hence the 
total current will be given by 


dl 

dz 


ffpre -<^> 2 

[2 -erf I 2 'l-erff L - z \l 

A 



dpdr 


(6.13) 


The first term may be integrated with the result 


/dl\ n - ea - 
\dz/j 2 /w 


(2ttR) 


(6.14) 


which is just the current per unit length for the case of no magnetic 
field. The other terms have been integrated numerically for a choice 
of values for R and p* characteristic of a Nike-Apache rocket in the 
ionosphere. The resulting currents normalized to the current of (6.14) 
are shown in Figure 15. 

It is apparent from the figure how for a typical Nike-Apache pay- 
load length of 60" the current per unit length at the mid-point of the 
rocket has decreased ty two orders of magnitude. In fact, the curve 
for L = 60" is essentially identical to the current par unit length due to 
electrons arriving from one direction only. This effect probably explains 
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NORMALIZED CURRENT PER UNIT LENGTH 



Figure 15. The Electron Current to a Cylinder in a Magnetic Field 
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why measured electron currents to probes on the nose-tips of Nike- 

Apache rockets have been observed to be low by a factor of two when 

95 

compared with normal diffusion theory. 

A similar procedure has been followed to compute the electron 
current to a sphere as a function of the ratio of the sphere's radius to 
the most probable radius of gyration, p*. The geometry is illustrated 
in Figure 16, where the magnetic field is parallel to the z-axis. The 
equations for the spival described by the electron are 


and 



(6.15) 


where (x t , ) are the co-ordinates of the guiding center. When this is 

combined with the equation for the sphere, r = R, the following equation 
is obtained for the value of z where the particle is incident on the sphere 


R 2 - z 2 = r 2 


+ 


P 2 + 2r t p 



(6.16) 


The left-hand and right-hand sides of this equation are also shown in 
Figure 16 as a function of z . It is apparent that all particles with 
( fj - p ) 2 > R 2 cannot be collected, while all particles with (i-j + p) < R 
will be collected. Those with (r t -p ) 2 < R 2 < (rj +p ) 2 may or may not 
be collected depending upon the particle's phase angle 4> and its wave- 
length. The critical phase angle, <f> c , is defined as the phase angle 
necessary so that the particle merely grazes the sphere. It is found 
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Figure 16. The Geometry of Electron Collection by a Sphere 

in a Magnetic Field 


72 



by combining Equation (6.16) with the equation obtained by differentiat- 
ing (6.16) with respect to z. We find 


4> c = arc cos 



2r x p 



CO z 


c 


V 

z 


(6.17) 


where 




(6.18) 


All particles with -0 c < 0 < + 0^ will not be collected. 

The case when 0 c = 0 occurs when the particle trajectory grazes 
the sphere twice, once when approaching and once when leaving. This 
defines a minimum wavelength, \ m = 2vm_ /Be (v ) , for a ^iven Tj and 

p, such that all particles with shorter wavelengths will be collected 
regardless of their phase angle. We obtain for the total current: 


I = 


4n_ ea_ 

ft P* 2 


r r r R - r i /*°° r 2?r 
" 0 % 'p=0 *'u=0 


Q(r r p, u, 0) d0dudpdfj 


+ 



Qd0du + 2 


00 


1 I Qd0du 

dpdr l > 

Ju=u m J<t> c 

J 


(6.19) 


where u = v z /a_, and where 


Q(r r p, u, 0) = 


upe - u 2 e -^/^) 2 


( 6 . 20 ) 


The area A is the same as shown in Figure 14. All the integrations 
may be carried out in a closed form except for the term containing 0 C . 
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The total current to the sphere in terms of the current with no magnetic 

field is 



(2-erf R/p*) ♦ <5^5 (l 


e - ( «w + »F(»/p , ) y 


( 6 . 21 ) 


where is the result of a computer integration of the term in 

(6.19) containing <f > c . This current is plotted in Figure 17, verifying the 
effect of a magnetic field in reducing the electron current to a large 
sphere by a factor of two. 
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Figure 17. The Electron Current to a Sphere in a Magnetic Field 
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CHAPTER VII 


EQUILIBRIUM POTENTIALS FOR VARIOUS 
ENVIRONMENTS 

1. General Equations and Computational Procedure. The general pro- 
cedure for finding the equilibrium potential of a body is to set the total 
current to the body equal to zero. The corresponding potential is the 
equilibrium potential, in accordance with the discussion in Chapter I. 

In this section some general results are presented in the form of equa- 
tions and graphs showing expected equilibrium potentials as a function 
of various parameters such as temperature, body size, body velocity, 
etc. In the following sections equilibrium potentials will be compu d 
for representative conditions in the upper atmosphere of the earth and 
in interplanetary space. 

For a small spherical body at rest we use (see Table 1) Equation 
(3.3) for ions and Equation (3.2) for electrons and find 

e ^/ kT - =1-^+0 (7.1) 

\ m - T +/ kT + 

as long as the body is negative. Here Q is the normalized current 
density 

0 = (7.?' 

2n e a + 

where J is the sum of both photoemission current density and current 
densities of energetic particle fluxes and their secondaries. This is a 
transcendental equation which must be solved numerically. Solution 
for Q = 0 are shown in Figure 18 as a function of the temperature ratio 
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NORMALIZED EQUILIBRIUM POTENTIAL,^ e/kT_ 



Figure 18 . Equilibrium Potentials for a Small Body at Rest without 
Photoemission or Secondary Emission 
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T„ /T + for ion environments of 0 + , He + and H + , As one would expect, a 

body is more negative when the atmospheric ion is heavier because of 

the lower thermal velocity. Although the normalized potential becomes 

less negative for an increasing temperature ratio, the potential itself 

becomes more negative. In other words, if T + were constant but T_ 

were increased, the potential of the body would become more negative 

even through the quantity <£e/kT. becomes smaller in magnitude. There 

are both theoretical and experimental reasons for believing that thermal 

96 

equilibrium is not always present in the atmosphere. 

Equilibrium potentials as a function of Q are shown in Figure 19 
for 0 + and H + environments and for three temperature ratios. At large 
values of Q the curves become quite flat. This is because the ion cur- 
rent is becoming unimportant in comparison with the (positive) photo- 
and/or energetic particle currents. When the ion current is negligible 
we find 


4>e „ 1 
= In 

kT_ 


f {v J 

,2n e a 


(7.3) 


as long as 



■ 4) 


When 0 exceeds the right-hand side of this inequality the body becomes 
positive. When this is true Q is at least /m^/nT = 42.87. Hence the 
positive ion current may be disregarded, and we find 


0e/kT„ = Q p f , (<t>) + Q + [1 + Y + f 2 (<j>)] + 0. [1 - Y, f 2 ($)) - 1 (7.5) 
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H+T./T 
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Figure 19. Equilibrium Potentials for a Small Body at Rest as a 
Function of Photoemission and/or Secondary Emission 



where Q p , Q + , Q_ are the normalized photoe.mission and energetic par- 
ticle current densities. Y ± is the effective secondary electron yield for 
the energetic particles; and f^f^and f 2 (£) are curves of current versus 
voltage obtained for a small body from Figures 7 and 11 respectively, 
and normalized to unity at 4> - 0. 

Although a "small body at rest" may not at first glance seem to 
apply to many real situations, this is not the case. Micrometeorites or 
other dust particles have been shown to be important in and below the 
D-region of the ionosphere in providing a sink for atmospheric ions and 
electrons. Their potential which must be known to estimate the magni- 
tude of this effect may be calculated from the preceeding equations and 
graphs. The equilibrium potential of small particles in interplanetary 
space where the positive currents to particles are due to photoemission 
and solar wind protons may also be computed from Equations (7.3) or 
(7.5). 

For larger bodies sheath effects must be taken into account. As 
long as 0 = 0 the equilibrium potential of the body will be negative. 

The appropriate equations are (3.2), (3.4), (3.5) and (3.6). Solutions for 
the equilibrium potential have been obtained numerically and are shown 
in Figure 20 as a function of the ratio of the body radius to the Debye 
length. Magnetic effects have not been included. When the ratio r/L 
is large. Equation (3.7) may be used with the result: 


<pe 

kT. 


= - 1/2 In 



(7.6) 


The magnetic effect of Chapter VI, Section 2, may easily be included in 
this equation by subtracting from the right hand side the quantity In (I /I 0 ) 
of Figure 17. 
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- 2.0 
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Figure 21 shows solutions for the equilibrium potential of a small 
moving body as a function of its Mach number, M, for the case when 
Q = 0. The appropriate equations are (3.14) for ions and 3.2 for elec- 
trons as long as the potential is negative. We find 




<t>e/ kT_ v/77 

■' =t m 


(erf M) M 2 - 



(7.7) 


As the velocity initially increases the ion current first decreases, as 
ions cannot catch up to the satellite from behind, and then increases as 
more and more are swept up in front. As a result, the equilibrium 
potential first becomes more negative, reaches a maximum, and there- 
after becomes more positive. 

The effects of photoemission on a small moving body are shown in 
Figure 22. Note the change from a linear to a logarithmic scale for 
the abscissa at Q = 10. Equation (7.7) was used for negative potentials 
and the following equation (from (3.9) and (3.3)) for positive potentials 
when Q is large: 


1 

f,(0) 



4 


} 


= Q 


( 7 . 3 ' 


where G(M, <j>) is the curly bracket of Equation (3.9). The equation has 
been written in this form because it has been found that the easiest 
computational procedure is to compute Q for a series of choices for 

When 4> is positive the Mach number becomes relatively unimpor- 
tant. This again is because the ion current is small compared with the 
large positive photocurrent. The equilibrium potentials are almost 
independent of the nature of the ion. This can be verified by noting that 
the curves for 0 are displaced to the right by ? factor of four from 
those for H* , which is the ratio of the corresponding thermal velocities 
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NORMALIZED POTENTIAL, <f>e/kj. 
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Figure 21. Equilibrium Potentials for a Small Moving Body 

without Photoemission 
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Function of Photoemission 



and hence also the ratio of the corresponding definitions for Q. Since 
the electron current to a small body is unaffected by the body's motion, 
the equation for large 0 reduces to (7.5). The reason that the plasma 
temperature must appear explicitly is that the function giving the 

relative number of photoelectrons that escape is independent of T. 

Turning now to the case of a large body in motion, we must take 
into account the effect of the wake on the electron current as discussed 
in Section 3.4. Figure 23 shows computed equilibrium potentials as a 
function of Mach number for two bodies whose radii are 5 and 10 Debye 
lengths. In this figure Q is zero and no magnetic effects have been in- 
cluded. The ion current was computed according to the method described 
in the appendix. In addition, results are shown when the ion current was 
computed from the sheath approximation of Equation (A13) together 
with Equation (3.13). It is apparent that the two methods agree very 
well for Mach numbers greater than about three. 

Figure 24 shows the effect of photoemission in driving a large body 
to positive potentials. The radius has been taken to be 33.3 Debye lengths. 
That the equilibrium potential does not depend strongly on size once r/L 
exceeds 10 can be seen by comparing the potentials at Q = 0 with the 
corresponding solutions in Figure 23. The effect of the induced poten- 
tial gradient due to the magnetic field in driving the body more negative 
is shown by the dashed curve for a Mach number of unity and an oxygen 
ion atmosphere. In Figure 24 the sheath approximation for the ion cur- 
rent of Equation (A13) in the appendix was used with Equation (3.13) for 
negative potentials and Eq. (3.9) for positive potentials. 

2. Expected Equilibrium Potentials in the Ionosphere for Large and 
Small Bodies. To illustrate the range of equilibrium potentials expected 
in the earth's ionosphere, two model atmospheres have been chosen 
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Figure 23. The Effect of size on the Equilibrium Potentials of a Moving Body 
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as a Function of Photoemission 



corresponding to low and medium temperatures. The data in Chapter 

VIII will describe equilibrium potentials for high temperatures. Figure 

25 shows the ion densities and temperatures chosen for the cold model, 

and Figure 26 the corresponding quantities for the warm model. The 

procedure in constructing each model was as follows: A temperature 

for the isothermal region above 500 km was first chosen. The peak 

electron density was then chosen to be 2 X 10 6 cm" 3 at 280 km for the 

warm model and 2 X 10 s at 350 km for the cold model, in accordance 

96 

with typical diurnal maximum and minimum densities. Chandra's ex- 
pressions for electron density in the upper F region were then used to 

97 

compute densities from 200 km to 425 km. Below 500 km T„ was 

set equal to T + for the cold model and normalized to the neutral gas 

temperature variation as published by the US Committee on Extension 

98 

to the Standard Atmosphere. For the hot model the ion temperature 
was set equal to the standard atmosphere neutral gas temperature, but 
the electron temperature measured on NASA rocket 6.04 was used. This 
rocket experiment flown on March 26, 1961, found an electron temp- 
erature at 360 km of 1600°K. The ion densities above 400 km geo- 
potential altitude (425 km true altitude) were computed from Bauer's 

ternary ionosphere model, which uses that altitude as a reference 

, . 100,101 

level. 

Bauer's model was used up to 6000 km altitude. Temperatures irk 

the cold model were kept at 700° all the way out to 20,000 km altitude 

(4.14 earth radii), whereas for the warm model the temperature was 

102 

increased in accordance with recent measurements by Serbu. Serbu's 
observation that the electron density beyond 2.2 earth radii may be ap- 
proximated by a power law with an exponent of -3.4 was used to com- 
pute the total electron- density out to 20,000 km. 
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Figure 25. Cold Model Ionosphere 
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Figure 26. Warm Model Ionosphere 
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TEMPERATURE (ev) 




The equilibrium potential of both a small and a large body have 
been computed for both model atmospheres under conditions of darkness 
and sunlight. The body was assumed to have a velocity equal to that 
necessary for a circular orbit at the corresponding altitude. The mag- 
netic induction effect was included for the large body by assuming that 
it was moving perpendicularly to the earth's magnetic field. This would 
be true for a body moving in the earth's geomagnetic equatorial plane. 
Finally, the secondary emission yields for the three kinds of ions en- 
countered as well as for electrons (incident on Al) were included in the 
calculation. 

The equations used in the calculations will not be given here. They 
are similar to those given in the preceding section but more complicated 
in that the ion current in general involves a summation over the three 
ion constituents. The photoemission current density was taken to be 
4 X 10" 9 amp/cm 2 , and the radius of the la^^e body to be 1.0 meters. 

The results of the calculations are shown in Figures 27 and 28. 
Several general observations may be made. First, the potentials for 
the large and the small bodies are very similar, lying in the range from 
-0.1 to -0.6 volts, until the effect of photoemission begins to predominate. 
Photoemission begins to have a significant effect at about 1000 km al- 
titude, but positive potentials are not obtained until about 5000 km in a. 
cold atmosphere and 10,000 km or above in a warm atmosphere. In 
darkness the satellite potential is mainly a function of the temperature 
ilone, as can be seen from the constancy of each solid curve in the 
isothermal region of the atmosphere. At higher altitudes — above 10,000 
km - there u a much greater range in the possible values for the 
potential. 
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Figure 27. Equilibrium Potentials for a Small Body 
in Cold and Warm Model Ionospheres 



EQUILIBRIUM POTENTIAL (volts) 



Figure 28* Equilibrium Potentials for a Large Body 
in Cold and Warm Model Ionospheres 
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The effect on the equilibrium potential of energetic particles trapped 

in the earth's magnetic field has not been investigated in detail, except 

103 

for some estimates made recently by Kurt and Moroz. Figure 29 
exhibits typical energetic particle fluxes in the earth's magnetic equa- 
torial plane. Not shown is the inner zone proton belt between approxi- 
mately 1000 and 10,000 km consisting of 10 Mev and higher energy par- 
104 

tides. The maximum omnidirectional flux of these protons is about 

5 X 10 4 cm" 2 st ' -1 at about 4000 km, which is much too small a flux to 

be a significant charging current. Much more significant is the belt of 

105 

low energy (5KeV) protons found by Freeman and also described by 
Hilton et al.*^ Extremely large fluxes were observed; but the spacial 
extent and altitude of maximum flux of these protons is not known. These 
protons would have a significant effect on the equilibrium potential, as 
is shown by the curves marked with an "E" between 1100 and 1550 km. 
The effectiveness of these protons is enhanced by the production of 
secondary electrons with a yield of 2.8, corresponding to an aluminum 
surface. 

In the same region of space as the high energy proton belt there 
are known to be large fluxes of energetic electrons. There are appar- 
ently large fluctuations in these fluxes with time, but the variation and 

extent of these fluxes are not well known. A typical observation of elec- 
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trons with energies above 20 KeV is shown. At these energies the 
total secondary emission coefficient is 0.3 or less on aluminum. The 
effect of these fluxes on the equilibrium potentials shown in Figures 27 
and 28 was computed and found to be either insignificant or just barely 
distinguishable. 
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More significant are the outer zone protons with energies above 
100 KeV described by Davis and Willianison. They found an expon- 
ential energy spectrum with a scale factor of 400 KeV at the altitudes 
shown* An integration of the secondary electron yields of Figure 13 
over this energy spectrum results in an effective yield of 2.25 for these 
protons incident on aluminum, and 1.34 for incidence on tungsten. The 
effect of these protons on the equilibrium potential is shown above 
10,000 km for the case of an aluminum surface in darkness. A body 
with a tungsten surface would have a potential between the two curves. 

In the sunlight the effective "Q + " due to these protons is at most less 
than 4% of the value of Q due to photoemission. 

It is emphasized that these choices of fluxes of energetic particles 
are illustrative only. In reality the fluxes vary widely not only in time 
but also as one moves away from the magnetic equator. This discussion 
is primarily for the purpose of demonstrating that the energetic particles 
trapped in the earth's magnetic field can have a significant effect on the 
equilibrium potential of a body. 

3. Expected Equilibrium Potentials in the Earth's Magnetosphere and 
in Interplanetary Space. The earth's magnetosphere may be defined as 
that region of space where the motion of charged particles is controlled 
by the earth's magnetic field. It is characterized by a relatively hot 
ionized gas with a temperature from a few thousand to 50,000°K, and by 
large fluxes of energetic electrons. ^®9»110 -j'jjgj.g j s a we u defined 
outer boundary to the magnetosphere which is quite analogous to the 
hydrodynamic description of the supersonic flow at a fluid around a blunt 
object. In this case the blunt object is the magnetic field and the fluid 
if, the expanding ion and electron gas coming from the solar corona, 
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aptly called the solar wind. At the magnetosphere boundary (magneto- 
pau e) there is a sharp drop in the magnetic field, a narrow transition 
region, and then a shock-front caused by the solar wind impact. Inter- 
planetary space is characterized by the solar wind streaming radially 
outwards from the sun at velocities from 300 to 10 3 km/sec, and with 
densities of 5 to 20 ions /cm 3 1 ^ 

The calculation of equilibrium potentials in both these regions is 

£ 

considerably simplified by the fact that the effect of low energy ions may 
be usually disregarded. Within the magnetosphere the flux of energetic 
(> 30eV) electrons is apparently usually larger than the total positive 
ion flux — at least at distances beyond about six earth radii from the 
earth's center. Consequently, the equilibrium potential of a body in the 
magnetosphere will be determined by a balance between the effects of 
photoemission and incident electrons only. 

Another simplification is the fact that the low particle densities 
and high temperatures result in a large value for the Debye length — 
typically 2 meters or larger. Therefore, it is not a bad approximation 
to use the small body equations for bodies as large, even, as typical 
spacecraft. Finally, because of both the negligible effect of low energy 
ions and the relatively low velocities of bodies in these regions with 
respect to the electron thermal velocity, the body may be regarded to 
be at rest. 

The equilibrium potential in the magnetosphere is very sensitive 
to the ratio of energetic electron flux to photoemission flux. For a 
negative body we find 


kT_ 

In 

e 



1 - 



(1 



(7.9) 
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as long as 


J.. J_ 2n ea. 

y (1 -y.) <1 <y(l -YJ + 

j p J P \G J p 

When the energetic electron flux dominates photoemission the potential 
may rise to a very high negative value — to a potential characteristic of 
the energetic electrons* In this case the low energy positive ions should 
be included in the calculation. However, this is not a very likely possi- 
bility for bodies with metallic surfaces. Maximum observed energetic 
electron fluxes are on the order of 10 9 /cm 2 -sec, yielding a current of 
~ 10” 10 amp /cm 2 , which is smaller than expected photoemission cur- 
rents. In addition, these large fluxes were for electrons in the 
200 eV to 40 KeV energy range, so that secondary emission is likely to 
be significant in reducing the efficiency of those electrons. 

When the right hand side of this inequality is violated the body will 
become positive: 


0e _ J P 
kT_ 2n e a. 



(7.10) 


Solutions to (7.10) are shown in Figure 30 for an electron energy spec- 
trum with two peaks: One at an energy (~500 eV) such that the second- 
ary yield is unity, and one at a higher energy where the yield is 0.3. 

The potential is limited to a few volts positive as the photo- and sec- 
ondary electrons are returned to the body. It was found that when the 
ratio J 0 /J p » corresponding to the higher energy flux, exceeded 0.14, 
the potential did not exceed +1 volt for the values of J i /J p shown in- the 
figure. 
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Figure 30. Positive Equilibrium Potentials for a Body in the Magnetosphere 




In interplanetary space the protons in the solar wind should be taken 
into account, although they contribute only about 10“ 10 amp/cm 2 , a factor 
of 10 less than photoemission. However, their effectiveness is not re- 
duced much at positive potentials as is the case for photoemission. 
Electrons in the solar wind apparently have two distinct energies — one 

from 3-5eV with densities equal to the proton density, and a higher tem- 

113 

perature group (20-50eV) with an order of magnitude smaller density. 
The equilibrium negative potential is 


■*• = — In [J, (1 + Y,) + J » (7.4) 

e . 1 2n e a_ p J 

The low temperature values for T_, n_ , u_ shox Id be used as long as 
the potential is not so negative that the flux of low temperature electrons 
has been reduced to a value less than that of the more energetic group. 
When this does occur the values for T_ , n_, a_ should be that of the 
more energetic electrons. 

In general one would expect a positive equilibrium potential in inter- 
planetary space unless the material of the body is such that photoemis- 
sion is unimportant. Solutions of the following equation: 



kT_ 2n_ e a_ 


f J + 

c Jp 


[i +v + f 2 (0m -i 


(7.12) 


are shown in Figure 31 for values of potential between 0 and +20 volts. 
The secondary yield has been neglected for proton energies less than 
1 KeV, and has been taken to be 2.0 for the 3KeV protons. The latter 
energy is expected only occasionally, as after large solar flares. Typ- 
ical values for Q_ and J + /J p are 10 to 100 and 10” 1 to 1 respectively. 
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Figure 31. Positive Equilibrium Potentials for a Body in Interplanetary Space 



It is apparent that a considerable positive potential cculd develop. Once 
+?.0 volts is exceeded the secondary electrons cannot escape, and the 
equilibrium potential is simply 



kT_ 2n e a_ 


(7.13) 
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CHAPTER VIII 


EQUILIBRIUM POTENTIALS MEASURED ON THE 
EXPLORER VIII SATELLITE 

1* Description of Experiment. The Explorer VIII Satellite was launched 
on November 3, i960, from Cape Canaveral, Florida, into an orbit with 
an inclination to the equator of 50°, a perigee of 420 km and an apogee 
of 2300 km. Its primary mission was the direct measurements of elec- 
tron density and temperature, positive ion concentration and mass, and 
the interaction between the vehicle and the ionized atmosphere. To ac- 
complish this mission the following experiments were flown: three cur- 
rent monitors were situated on the satellite equator, two consisting of 
single-grid charged particle traps appropriately biased to measure 
positive ion and electron current densities, and one consisting of an ex- 
posed plate to monitor the total net current. Two electron traps to 
measure electron density and temperature were mounted on the top cone 
of the satellite (see Figure 32); and one ion trap to measure positive 
ion density, temperature and mass was also mounted on the equator. In 
addition, there was an electric field meter mounted on the top surface 
and a radio -frequency plasma impedance probe consisting of two 10 foot 
wires extending outward from the equator along opposite radii. In addi- 
tion to the preceding experiments, there were also two micrometeorite 
experiments and an attitude sensing system consisting of a solar cell 
and a horizon sensor, both mounted on the equator. 

The data from the experiments were sent by a 108 Me telemetry 
system to any of the eleven NASA ground stations within range when 
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the data were obtained. All the instrumentation operated continously 
except the electric field meter which was turned on for two minute 
periods by radio command from the ground. 

The shell of the satellite was made of aluminum and consisted of 
two truncated cones joined at the equator by a short cylinder 30 inches 
in diameter and 6-1/8 inches high. The axis of the satellite was also 
30 inches long, including the 3-inch cylinder at the bottom used for con- 
nection to the booster. The total exposed surface area was 2750 in 2 , 
which yields an effective radius for an equivalent sphere of 14.8 inches. 

Shortly after launch the satellite was "de-spun" to a spin-rate of 
approximately 20 rpm. The celestial co-ordinates of the positive spin- 
axis remained in the region of zero degrees declination and 180° right 
ascension during the satellite's useful life of five weeks. The data from 
the horizon sensor proved to be difficult to analyze, with the result that 
there was an uncertainty of about ±15° in the spin-axis co-ordinates. 

The angle between the spin axis and the sun was known to remain close 
to 60°. 

A more detailed description of the satellite and the instrumentation 
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is available elsewhere. Some of the results of the experiments which 

have been reported are measurements of the sheath currents, including 

the first experimental verification of the induced potential gradient due 

40 

to the satellite's motion in the earth's magnetic field; the first direct 
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detection of helium ions in the earth's upper atmosphere; and the 
diurnal variation of temperature in the upper atmosphere. 

This investigation is concerned only with the results obtained from 
the si gle-grid electron trap. A block diagram of the experiment is 
shown in Figure 33. The sensor consists of a grid flush with the satellite 
skin with a sweep potential varying from -1.2 to +8 volts and back in 0.4 


107 










sec as shown in the figure inset. Behind this grid is a collector biased 
at +15 volts to remove photoemission and incoming ion current from 
the measured collector current. The electrometer employed 100% neg- 
ative feedback to maintain stability against drift and to keep the collector 
at a constant potential independent of the measured current. It was time- 
shared with the ion current monitor and the ion trap, each experiment 
being connected for a 30 second period consisting of three 10-second 
intervals on each of three electrometer ranges. Currents from 10“ 9 
amp to 5 X 10" 5 amp could be measured. 

The principle of the experiment is as follows: as long as the grid 
is negative with respect to the plasma, the current decreases exponen- 
tially with voltage in accordance with Equation (3.2). The slope of the 
straight line obtained when the logarithm of the current is plotted against 
potential is determined by the electron temperature: 


d(ln I) e 
dV " kT 


( 8 . 1 ) 


When the grid is positive with respect to the plasma, electrons are 
attracted and the current increases at a slower rate which depends on 
the geometry as well as the electron temperature. The potential V s at 
which this change in character of the electron current occurs identifies 
the time when the grid was at the same potential as the plasma. Since . 
the measured potential is with respect to the satellite, we identify 4> K > 
the satellite potential with respect to the plasma, as - V s . The electron 
density may be simply computed from the current at that point by the 
relation 

/kT v /a 

(I) v.v. ■»•** (sSj <8 ' 2 ' 


109 



where t is the electrical transparency of the grid (53%), and A is the 
aperture area of 13,0 cm . 

Theoretically it should also be possible to obtain the satellite poten- 
tial from the shape of the current versus voltage curve obtained from 
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the ion trap. Actually, this proved to be quite difficult because of 
the strong dependence of the current on the angle between the normal to 
the trap and the satellite velocity vector. 

2. Experimental Results . Figure 34 shows a typical current-voltage 
curve obtained with the electron trap. The circled points were obtained 
with the electrometer in the medium sensitivity range and the crosses 
with the high sensitivity range. In this case the data from the two ranges 
was obtained 1.5 seconds or one half of a spin period apart. In general, 
the data was chosen to be either a full period apart or very close to- 
gether on either side of the range switching, although this was not always 
possible. Successive points are separated in time by about 0.010 sec- 
onds. The currents below about 2 X 10“ 8 amp were corrected for a 
displacement current effect due to the changing voltage on the grid. 

The charr - teristically linear portion of the curve on this semi-log 
plot with a distinct change in slope at 5.2 X 10~ 7 amp is plainly apparent. 
The following procedure was used to compute the temperature, density 
and satellite potential: All currents greater than 1% of full scale on the 
low and medium sensitivity ranges and 10% of full scale on the high 
sensitivity ra*v*, but less than the apparent break-point in the slope, 
were fitted to a straight line by a least-squares calculation. This line 
was plotted and the fit to the data was examined at and above the break- 
point to see if points there should be added or subtracted. Then a new 
least-squares calculation was made with the new set of points. This 
process was repeated until there was no question that the best possible 
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Figure 34. Typical Current- Voltage Curve Obtained 
with the Electron Trap 



estimate of slope and break-point had been obtained. The density, sat- 
ellite potential and electron temperature with its standard deviation 
were then computed ham Equations (8.1) and (8.2). 

This procedure was used to compute densities, temperatures and 
satellite potentials from some 250 current-voltage curves. The tem- 
peratures and densities obtained on magnetically quiet days (according 
to the Committee on Characterization of Magnetic Disturbances**®’***^) 
are presented in Figures 35 and 36. The satellite potentials are pre- 
sented in Figures 37 through 39. The data has been separated according 
to magnetic activity because there was a significant difference in the 
results between magnetically quiet, medium or disturbed conditions. 
During a given condition, however, there was no significant difference 
between results on different days that was not obscured by the experi- 
mental scatter. 

It is reasonable to expect atmospheric conditions to stay reasonably 
stable from day to day over the Explorer VIII orbit for a given condition, 
because of the fact that the local mean time was practically a fixed func- 
tion of position in the orbit. Thus, diurnal variations would not appear 
explicitly but would be folded into the variation with altitude in the same 
way from orbit to orbit. This feature of the or j;.t is of course due to 
the fact that the plane of the orbit is fixed in inertial space and conse- . 
quently rotates quite slowly (one degree per day) with respect to the 
sun. The local mean time shown in Figures 35-39 is for November 20, 
and is good to within about ±1 hour over the active life of approximately 
one month. 

For these reasons, plus a consideration of the sources of error 
in the measurements, if is felt that the scatter in the data is 
primarily experimental rather than a reflection of real geophysical 
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igure 35. Measured Electron Temperatures for Magnetically Quiet 
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Figure 36. Measured Electron Densities for Magnetically Quiet Days; 
Line Indicates Values Used in Calculations 
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Figure 37. Measured and Calculated Equilibrium Potentials for 
Explorer VIII on Magnetically Quiet Days 
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Figure 38. Measured and Calculated Equilibrium Potentials for 
Explorer VIII on Moderately Magnetically Active Days 
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Figure 39* Measured and Calculated Equilibrium Potentials for 
Explorer VIII on Magnetically Disturbed Days 
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variations. The typical standard deviation for the measured tem- 
perature was about 25% rather than the 11% of Figure 34. The 
measured satellite potential is in turn very sensitive to the meas- 
ured slope of the current-voltage curve which determines the position 
of the break. The absolute value of the grid voltage is known within an 
uncertainty of about ±0.08 volts. This uncertainty affects the determin- 
ation of $ s more strongly than it affects the determination of the tem- 
perature where only the relative voltage is needed. 

The uncertainty in the measured currents is about 1% of full scale 
for the low and medium sensitivity ranges and 10 % for the high sensitiv- 
ity ranges, as indicated previously, c computed density reflects not 
only this uncertainty but also the uncertainties in the temperature and 
in the determination of the break-point. 

Before the measured satellite potentials can be compared with 
theoretical values, the dimensionless potential 7 ? s = <£ a e/kT must be 
computed. The relative error in t) % is the sum of the relative errors 
in and T . Consequently, it was decided that a useful comparison of 
the measured satellite potentials with theoretical values could best be 
done by using a smooth curve of temperature and density over the orbit 
to compute the expected potentials. This way the non-systematic dis- 
crepancies between the measured and predicted values reflect primarily 
the scatter in the measured 4> a alone. The smoothed temperatures 
representing a running mean for the three magnetic conditions are shown 
in Figure 35. These were used both in the theoretical calculation of 17 
and in the normalization of the measured . Only the electron density' 
data for the magnetically quiet days is shown in Figure 36. 
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The following expression was used to compute the expected values 


for 7] : 


(i-s) m; 
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A.V^T Jp 
Ay2n ea + 


(8.3) 


In this expression S and Y are the reflection coefficient and sec- 
ondary yield for electrons and oxygen ions respectively. Bauer has 
shown that for the high temperatures found here, oxygen ions will be 
predominant over the whole orbit. L is the Debye length and R the 
effective radius of the satellite of 0.37 meters. The quantity (L/R) 2 is 
used in accordance with the discussion in Section 3.4 on the effect of the 
wake. The term containing p represents the magnetic induction effect 
of Equation (6.3), and F(R/p*) is the ordinate of Figure 17. The cross 
product of V and B is a fixed function of position in the orbit if it is 
assumed that the earth's magnetic field is a dipole with an axis along 
the earth's geographic axis. M is the Mach number, V/a + ; t ( 77 ) the 
sheath thickness after Equation (A13) of the appendix; and Ay , A s , A y , 
are the total area and the projected areas in the direction of the sun 
and velocity vector respectively. It has been assumed that the ion tem- 
perature is equal to the electron temperature. 

The results of the computations are shown by the lines in Figures 
37 through 39* The dashed curve is for a photoemission current density, 
Jp, of 4 X 10~ 9 amp/cm 3 and the solid curve for 8 X 10~ 9 amp/cm 2 . 
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There is reasonable agreement at the higher altitudes between the curves 
and the data points, particularly with the higher photoemission current 
value. A systematic overestimate of the electron density has the same 
effect, on the predicted potential as an underestimate of photoemission. 
This is a more likely explanation of the better fit of the solid than the 
dashed line. For the disturbed days an even higher values for (J p /n) is 
indicated. 

Below 700 km the measured values of 17 are much more negative 
than is expected. The change in potential due to passage from sunlight 
to darkness is reflected by the change in the curves between 900 and 
1100 km on ascent and 900 and 700 km on descent, and is not of sufficient 
magnitude to be an explanation. 

It is suggested that the more negative potentials at low altitudes 
are caused by the rf impedance probe experiment on the satellite. This 
probe, consisting of the two 10-foot wires described previously, was 
operated continuously at a frequency of 6.5 Me. The amplitude of the 
rf voltage was about 0.3 volts. At high altitudes where the electron 
density is low, the local plasma frequency is much lower than 6.5 Me. 
Consequently, the additional current to the probe due to the rf is negli- 
gible, as is discussed in Section 5.7. As the satellite descends in altitude, 
the plasma frequency increases; 



and as the satellite approaches perigee, the plasma frequency will in' 
general go through the probe frequency. For example, the electron 
density of 10 6 /cm 3 occurring at 425 km in Figure 36 corresponds to a 
plasma frequency of 9 Me. Therefore, a very large increase in electron 
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current to the probe results as the probe frequency goes through the 
resonant frequency, which drives the whole satellite to a more negative 
potential. 

To put this suggestion on a quantitative basis, it is necessary to 
know the value of the current at the resonant peak and the resonant 
frequency, which is generally somewhat lower than the plasma frequency. 
Whale has shown that resonance occurs for a cylindrical dipole antenna 
at f /V5T, which yields a resonant frequency of 6.34 Me for the density of 
10 6 electrons/cm 3 at 425 km.^^ At these altitudes the resonant peak 
height is determined by a phase -mixing damping rather than collisional 
damping. 

This mechanism has not yet been studied extensively, but Crawford 

has suggested that the ratio (Ai /i Q ) peak to (Ai /i 0 ) f <<f p is probably 

121 

between 5 and 100 under these conditions in the ionosphere. Figure 

40 illustrates the additional current to a probe as a function of f / f p . 

94 

This curve was computed from Equation (19) of Harp and Crawford 
with an assumed collision frequency to give a peak current ratio of 20. 
The resonant frequency here is at 0.68 f . 

At 425 km for the quiet model, where n_= 1 X 10 6 /cm 3 , the in- 
creased current to the probe at f / f p = (6. 5/9*0) Me = 0.72 is Ai = 39 i Q , 
where i 0 is the current to the probe with no rf. When this added elec- 
tron current is taken into account in the left-hand side of (8.3), the new 
equilibrium potential, (0 # eAcT), is -4.16 rather than -2.1. This potential 
is indicated by the curve marked (I) in Figure 37. The value of Ai/i Q 
necessary at 500 km to give the indicated potential of -3.3 is 11.9 where 
f / f p is approximately 1.6. 

It should be emphasized again that the shape of the resonance curve 
.in Figure 40 is based on the collisional damping model, which does not 
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Figure 40. The Response of a Radio -Frequency Probe 


in a Plasma 



apply at these altitudes in the ionosphere. There is evidence from some 
resonance probes flown by Japanese experimenters that the resonance 
peak in the ionosphere is much broader. This is also indicated by 
the Explorer VIII potential data. In general, the potential begins to go 
more negative between 50C and 600 km where the electron density is 
near 10 s /cm 3 rather than 10® /cm 3 , and hence f/f 2.3. However, 

p 

the fact that this effect depends so strongly upon the electron density 
means that the smoothed model for the density in the computations will 
not reproduce the measured potentials closely. 

The curves marked (2) and (3) have been computed by assuming 
that at 425 km the probe frequency is at the resonant frequency. The 
ratios of (Ai /i Q ) to (Ai /i 0 ) 0 have been assumed to be 28.5 and 144 
respectively, giving values for (0 s e/kT) of -4.59 and -6.12. The corres- 
ponding values for the current ratio at 500 km are 13 and 27. 

These calculations show that the current resonance effect as the 
plisma frequency approaches the probe frequency can quantitatively 
account for the negative potentials observed near Explorer VIII perigee. 

It is necessary that the resonant peak be broader than what has been 
observed in laboratory work; but following the indications of the Japanese 
experiments, it is suggested that the resonant effect is broader at low 
neutral particle densities where collisional damping does not occur. 

3. Conclusions . Two major conclusions may be drawn from these 
results. The first is that there is general agreement between the pre- 
dicted values of equilibrium potential and the measured values. There- 
fore, the various mechanism of charge collection that have been dis- 
cussed have been evaluated correctly as far as their importance for 
Explorer VIII is concerned. Further, no important mechanisms have 
been omitted. 
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The second conclusion is that much more experimental work re- 
mains to be done. An order of magnitude improvement in the accuracy 
of satellite potential measurement is required. And to evaluate the 
several charge collection mechanisms, precise measurements of the 
environmental plasma properties are required simultaneously. If further 
experimental study of satellite potential is warranted, it should prefer- 
ably be a satellite designed primarily to that end. This way both the 
surface characteristics and the geometry of the satellite may be designed 
with their effect on equilibrium potential in mind. Any feasible experi- 
ments that will contribute to an understanding of the potential should be 
included, and they should be carefully designed so that any effects they 
might have on the potential themselves can be controlled. 

Several areas of further laboratory investigation are also indicated 
by this study. Much work needs to be done on photoemission yields of 
materials in the ultra-violet wavelengths. The yields of aluminum 
especially need to be verified for various surface conditions. Magnesium 
should be investigated, and also dielectric materials such as mylar that 
are often used on satellite surfaces. Work is also needed on the sec- 
ondary yields for ion impacts, especially in the energy range from 1 to 
10 KeV. 

The theory presented in Chapter VI on the effect of a magnetic field 
on the collection of electrons is a subject that could easily be pursued 
in the laboratory, although the extension to the case where electric fields 
are present will undoubtedly be very difficult. 

Finally, the radio -frequency resonance phenomenon at low pressures 
nev i'i study both experimentally and theoretically. It may be that the 
best place to pursue this experimentally would be from a satellite where 
wall-effects can be avoided, and where probes large compared to the 
Debye length can be employed. 
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APPENDIX A 


ION CURRENT TO A NEGATIVELY -CHARGED 
MOVING SPHERE 


Poisson's equation in spherical coordinates with the space charge 
given by Equation (3.15) may be written 


+ l ^ = 1 _ e-v 
dx 3 x dx 


(Al) 


where y is the normalized potential. 



( A2) 


and x the radial distance normalized to the Debye length. 

Solutions to this equation subject to the boundary condition that 
y 0 as x -* oo are plotted in Figure Al. It should be noted that the 
solution corresponding to any interior boundary condition defined by the 
pair of values (x Q , y Q ), the potential o the sphere of radius r Q , is des- 
cribed external to x 0 by the single curve on which x Q , y 0 lie. 

The approximation made in assuming Equation (3.15) for the space 
charge is that the ion density is undisturbed, and only the electrons 
respond to the field in accordance with the Boltzmann factor. Thus, 
this equation applies to the case where the body's speed is much larger 
than the ion thermal velocity. In the satellite coordinate system the 
ions are approaching with uniform speed from one direction. All those 

with impact parameters less than a certain value will be collected. 

38 

Walker in his thesis has shown that is is possible to compute this 
impact parameter from consideration of the conservation of energy and 
angular momentum without actually obtaining the particle trajectories. 
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NORMALIZED POTENTIAL, y = -<£e/kT 



Figure Al. The Potential in Front of a Negatively-Charged Moving 
Sphere; Contours of the Minimum Distance of Closest Approach are 
also Shown for Various Ion Energies 
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Following his discussion, but generalizing to any functional form 
for the potential (as long as it decreases monotonically), rather than 
restricting ourselves to a power law as he does, we define the pitch 
angle a as the angle between the particle velocity vector and the radius 
vector. The pitch angle at any point for a particle with impact param- 
eter b and initial velocity V is given by 



(A3) 


where 0 is the value of the potential at the point r . 

We divide all trajectories into two classes: periastron trajectories 
are followed by particles which would have (if there were no absorbing 
surface) a distance of closest approach, r min , corresponding to the value 
of r where a = tt/2, Pericritical trajectories are orbits which spiral 
in towards the origin with a never reaching W 2. Such trajectories can 
only exist if the potential falls off more strongly than r** 2 . This can be 
seen from Equation (A3) where sin a always increases as r decreases 
along a trajectory unless 0 falls off faster than r” 2 . In this case sin a 
may have a maximum value less than 7r/2. 

For a given initial kinetic energy, ions with large b will always 
describe periastron orbits. If pericritical orbits exist then there will 
be some minimum impact parameter for which a may equal 7t/2. Al- 
ternately, there is some minimum distance of closest approach, * £t 
which a particle may arrive and still escape. To find this we differentiate 


b = r 


mi 


■f- 


2 < *< r »in> e 


mV 2 


(A4) 
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with respect to r Bin , and find that r c is defined by the point where 


+ § (&) ■ + & < A5 » 

Since 1/2 mV 2 = 1/2 im^ + ^e, this equation describes the surface where 
the centripetal acceleration is equal to the centrifugal force per unit 
mass, and it is clear why particles that cross the surface r c will be 
accelerated towards the origin. 

The corresponding critical impact parameter, b , which divides 
periastron from pericritical trajectories is given by 


or 


/ 2</>(r c )e _ A c e / d0 \ 

>c ~ r<J f " mV 2 r * f mV 2 (dr) 



(A6) 


in terms of normalized quantities, where u = mV 2 /2kT. 

The program which computed solutions to Equation (Al) also tab- 
ulated along with each solution the value of u for a series of values of 
x c in accordance with Equation (A5). The corresponding values of x bc 
were also tabulated. Some of these results are also shown in Figure 
Al in the form of contours of the distance x e for a constant initial 
energy u. It is apparent that in general each solution to (Al) has two 
points where Equation (A5) is satisfied. What this means physically can 
be seen from Figures A2 and A3. In these figures the equivalent one- 
dimensional potential, Y(x) , is plotted against x , where Y(x) is given by 

Y(x) = <x b /x) 2 u - y(x) 
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(A7) 




Figure A2. Equivalent Potentials for a Given Potential Distribution and 
Energy u, but for Various Impact Parameters 
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Figure A3. Equivalent Potentials for Various Potential Distributions 
Denoted by (x o ,y o )» but for a Given Energy and Impact Parameter 
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(See, for example, Goldstein's discussion of the "fictitious potential" 
for a two-body central force problem), 123 It is clear from the figures that 
true pericritical trajectories cannot exist because of the repulsive cen- 
trifugal barrier which always dominates the potential at small values 
of x. It is also generally true that the repulsive centrifugal barrier 
dominates at large values of x because of the charge on the body by the 
sheath. However, there may still be a potential well at intermediate 
values of x. If a particle is to penetrate into this region it must have 
at least an energy equal to the maximum height of the barrier outside 
the well. 

The distance of closest approach of a particle is the largest value 
of x in Figure A2 for which the equivalent potential crosses the hori- 
zontal straight line corresponding to the particle energy. It is clear, 
then, that the critical impact parameter, x fae , describes the case when 
the equivalent potential is just tangent to the particle energy line at the 
barrier maximum outside the well. The critical radius, x , which de- 

C 

scribes the turning point thus has two values. One is at the position of 
the barrier maximum and the others is the inner turning point — i.e., 
the boundaries of region B. Whether or not a particle is collected by 
the body depends on the relation of the body's radius to the position of 
the turning point. Particles with impact parameters greater than x bc , 
the critical impact parameter, will be reflected at the outer centrifugal 
barrier, while particles with impact parameters less than x bc will come 
on in to be reflected at the inner barrier. 

Three cases may be distinguished corresponding to which of the 
three regions A, B or C contains the body's radius x 0 . If x 0 is in re- 
gion A, a grazing trajectory exists for an impact parameter x g less 
than x bc , and is given by Equation (A4), 
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Xg = X 0 /I + y 0 /u 


(A8) 


If x 0 is in region B, a grazing trajectory does not exist (except for 
trapped particles). All particles with impact parameters less than x bc 
will spiral in and be collected so that the impact parameter for collec- 
tion is x bc . In region C grazing orbits again exist and the impact param- 
eter for collection is again x g of Equation (A8). The rules for collection 
may be summarized as follows: If x Q < x bc < x g , use x bc to compute 
the effective cross-section; otherwise use x g . The resulting current- 
voltage curves for various combinations of x fl and u are shown in 
Figure 4 of Chapter III. 

These rules depend on the fact that the various equivalent potential 
curves for a given u and (x fl , y 0 )» hut for various values of x fa , do not 
intersect. This can easily be verified by forming the difference y 2 (x) -y 2 (x) 
= u/x 2 (x£j -x b2 ) > 0 for all x, where x bl > x b2< 

Figure A3 illustrates how the current to a body changes as the po- 
tential is increased. At small potentials corresponding to the uppermost 
curve the effective cross-section is 7 tx 2 , which increased linearly with 
the potential. As the potential increases a well forms which increases 
in width and depth, so that eventually x fl moves into and remains in 
region B. The current is then 

i = >&/*? (A9) 

l 0 

which increases more slowly with the potential. Again, it is necessary 
that the various curves for a given u and x fa but for various solutions 
to (Al) do not intersect. This may be verified by forming Yj(x) - Y 2 (x) = 
y 2 (x) - yj(x). Since the solutions y shown in Figure Al do not inter- 
sect this is true of the equivalent potentials as well. 
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Following Gringauz, we may define x bc as the sheath radius. For 
large bodies we may make use of this analysis to derive an analytic 
expression for the sheath radius which may be used in Equations (A9) 
or (3.13). 

For large bodies the second term on the left of Equation (Al) may 
be neglected, and the right-hand side may be linearized in accordance 
with the discussion of Section 2 of Chapter I. The solution is then simply 


— (x— *o ) 

y = y 0 e 


(A10) 


and the condition for x is 

C 

(l-x/2) =-.!!_ 
Vo 

To a first approximation x c - x 0 is given by 



(All) 


(A12) 


for 



and the sheath thickness t by 


Jt 

L 


= * bc “ * e 




( Al 3) 


The error in this approximation increases with increasing <p 0 and de- 
creasing x Q . However a comparison with the computer solution for the 
case where x Q = 23 and mV 2 /2kT = 5 indicated that the approximate 
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value for t of (A12) was in error by less than 1%, 10% and 25% for 
values for (-<£ 0 e/kT) of 1,0, 2.4, and 4.7 respectively. The error in 
the computed current will be still less, by a factor of about (2 1 /r 0 ). 
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APPENDIX B 


SUMMARY 

A body in the upper atmosphere or in space will acquire an electric 
charge through various mechanisms such as impacts by ions and elec- 
trons or electron emission. The charge on a body influences its rate of 
charging with the result that an equilibrium charge, or potential, is 
reached such that the net current to the body vanishes. Knowledge of 
the equilibrium potential of a body in space is needed to determine the 
motionof micrometeorites, the dragon earth satellites, and to assess the 
behavior of certain experiments on satellites. The various treatments 
of the problem in the literature have generally been restricted to a 
consideration of only a few charging mechanisms. Little data has been 
obtained, and in no case have the measurements been analyzed in terms 
of the expected potential with an evaluation of all the possible charging 
mechanisms. 

The collection of ions and electrons from the environmental plasma 
as a function of potential depends upon the body's size, shape, and ve- 
locity. The appropriate equations for ion and electron currents to a 
sphere are available in the literature for bodies small compared to a 
Debye length. For large bodies, an estimate of the influence of the 
plasma sheath is required to determine the current in an attractive 
field. Walker's estimate of sheath thickness is used for bodies with 
velocities small to the charged particle thermal velocity.* Poisson's 
equation has been solved numerically for high-velocity spheres, and- the 
ion current obtained by an analysis of the ion's distance of closest ap- 
proach. A negative space charge is formed in the wake behind such a 
body which reduces thfe electron current. 
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Photoemission is an important charging mechanism for bodies in 
sunlight. Measurements of photoelectric yields in the near and far 
ultraviolet for aluminum and tungsten are reviewed and compared with 
photocurrents measured above the atmosphere. A photocurrent density 
near 4 X 10" 9 amp/cm 2 is indicated for both materials. Secondary 
electron emission upon energetic ion or electron impact may also be 
an important charging mechanism, especially in the earth's radiation 
belts. Secondary yields taking into account reflection and primary 
back-diffusion are presented for electrons and for ions for energies 
up to 10 MeV. 

There are several other charging mechanisms that in general may 
be neglected except possibly in special circumstances. These include 
the effects of cosmic rays, radioactivity, thermionic and field emission, 
collisions with dust grains, and the influence of radio -frequency electric 
/ields. A magnetic field can affect the equilibrium potential of a body 
m two ways: the motion of the body through the field induces a potential 
gradient which results in an increased electron current; in addition, the 
restriction of the electron to a spiralling motion along the field line 
decreases the effective collection area of the body. Expressions de- 
scribing this latter effect are derived for both cylinders and spheres. 

These charging mechanisms are evaluated for conditions in the 
upper atmosphere and in interplanetary space, and are combined into 
expressions from which the equilibrium potential may be determined. 

In the ionosphere where electron densities are relatively high, the equi- 
librium potential is typically a few tenths of a volt negative. The im- 
portant mechanisms are environmental ion and electron collection. The 
energetic particles in the earth's radiation belts have a negligible effect 


136 



2 

except for the large fluxes of 5 KeV protons observed by Freeman and 
tho 100 KeV protons beyond 2 earth radio observed by Davis and William-* 

3 

son. At higher altitudes the equilibrium potential may become positive 
in the sunlight as photoemission predominates over positive ion collection. 

In the earth's magnetosphere where electron densities are low, large 
energetic particle fluxes may occur. The potential is sensitive to the 
ratio of electron flux to photoemission, and may vary widely.. Positive 
values, which are more likely, are limited to a few volts, but large neg- 
ative values are possible for large energetic electron fluxes. In inter- 
planetary space positive potentials due to both photoemission and the 
solar wind protons are expected, unless the material is such that photo- 
emission is unimportant. 

The equilibrium potential of the satellite Explorer VIII has been 
measured between the altitudes of 420 and 2300 km in both darkness 
and sunlight by means of an electron trap experiment. The results are 
compared with expected values calculated from the expressions derived 
earlier for equilibrium potential. The calculations include magnetic 
field effects and photoemission, rnd make use of simultaneous measure- 
ments of the plasma density and temperature. Although there is con- 
siderable scatter to the experimental potentials, there is general agree- 
ment at higher altitudes with perhaps a larger photocurrent of about 
8 X 10“ 9 amp/ cm 2 indicated. At low altitudes the measured potentials 
are more negative than anticipated. This is attributed to the effect of 
a radio -frequency plasma impedance experiment carried on the satellite. 

It is shown that near satellite perigee where the plasma frequency ap- 
proaches the probe frequency of 6.5 Me, the resonant increase of elec- 
tron current to the probe can quantitatively account for the more nega- 
tive potentials. 
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